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ИЗ ПРЕДИСЛОВИЯ АВТОРА 


Эта книга написана на основе курса лекций, прочитанных 
автором ocentio 1954 года в Калифорнийском техиологиче- 
было изложение 


ском институте. Основной целью курс 
различных методов асимитотического вычисления интегралов, 
содержащих большой параметр, а также методов решения 
дифференциальных уравнений с помощью  асимптотических 
разложений. Выбор столь широкой области вопросов привел 
к тому, что изложение приняло несколько эскизный характер. 

Глава | содержит краткое введение в общую теорию 
асимптотических разложений. Эта теория изложена лишь 
в пределах, необходимых Aid того, чтобы дать теоретическое 
обоснование главной части курса; эта скромная глава никоим 
образом не может заменить систематического и полного изло- 
жения рассматриваемого вопроса, даниого недавно ван дер 
Корпутом. 

В главе П изложены наиболее важные методы асимптоти- 
ческого вычисления функиий, заданных в виде определенных 
интегралов. Эта глава многим обязана прекрасной брошюр? 
Коисона, посвященной данному вопросу. Из-за недостатка 
времени мы дали лишь краткое изложение некоторых из MC- 


TOAOB асимптотиче ского вычисления (интегрирование по частям, 
метод Лапласа, метод перевала, метод стационарной фазы) 
и совсем не затронули вопросов, связанных с двойными и 
кратными интегралами. 

Последние две главы посвящены решению обыкновенных 
линейных дифференциальных уравнений. В главе ИП «большой» 
величиной, по которой ведется асимитотическое разложение, 
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является независимое переменное. Изложение ограничено 
рассмотрением дифференциальных уравнений второго порядка, 
для которых бесконечно удаленная точка является иррегу- 
лярной особой точкой первого ранга, а определяющее урав- 
нение имеет два различных корня. В главе ГУ «большой» ве- 
личиной является параметр в дифференциальном уравнении. 
Изложение ограничено рассмотрением = дифференциальных 
уравнений второго порядка в вещественной области, причем 
аргумент пробегает ограниченный замкнутый интервал. Даны 
оба приближения Лиувилля и различные их обобщения на 
случай, когда интервал содержит точку перехода. 


Декабрь 1955 г. | А. Эрдейи 


ВВЕДЕНИЕ 


Часто бывает так, что для вычисления некоторой вели- 
чины можно использовать расходящийся бесконечный ряд, 
причем сама величина в некотором смысле является суммой 
этого ряда. Типичная ситуация такова: некоторая функция 
разлагается в функциональный ряд, причем приближение, 
даваемое несколькими первыми членами ряда, тем лучше, 
чем ближе независимое переменное к некоторому предельному 
значению (таким значением часто является со). Во многих 
случаях члены ряда сначала быстро убывают (тем быстрее, 
чем ближе независимое переменное к предельному значению), 
но потом члены ряда вновь начинают возрастать. Такие ряды 
иногда называют полусходящимися (Стилтьес), а вычислители 
часто говорят о сходящихся в начале рядах (Эмде); однако 
в математической литературе, как правило, применяется назва- 
ние асимптотические ряды (Пуанкаре). Мы увидим ниже, что 
асимитотические ряды могут как сходиться, так и расхо- 
диться. 

Рассмотрим пример, исследованный впервые Эйлером (1754). 
Очевидно, что ряд 


S(x)==1— Иа Зы... = 1х" (1) 


расходится при всех x40. Однако при малых значениях 
x (скажем, х=107'?) члены ряда сначала весьма быстро 
убывают и можно вычислить приближенное численное значение 
5 (х). Возникает вопрос: какую функцию представляют при- 
ближенно эти численные значения? 

Введем, следуя Эйлеру, функцию 9 (х) = х5(х). Формально 
дифференцируя ряд (1), получим 


g (x) = 1! — 2k + 3x? ie et Lg 
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откуда имеем 
ху (хх) =x. 


Таким образом, %(х) можно истолковать как решение этого 
дифференциального уравйения, обращающееся в нуль при 
х=-0. С лругой стороны, используя эйлеров интеграл вто- 
рого рода 

ao 


ni=={ ee" dt, 


получаем 
nm 
e-'tdt4 as Ca eee = 


5 


$) = \е'4Ё—х 


® 8 
1-38 


= yer ryt ео, 


5 


Если просуммировать формально под знаком интеграла, то 
функция S (x) примет вид 


es 
{ dt. (2) 
о 


Но функция 


(3) 


является вполне определенной аналитической функцией от x 
в комплексной плоскости, разрезанной вдоль отрицательзой 
действительной полуоси, причем эта функция тесно связана 
с так называемой интегральной показательной функцией Ei (x). 
Возникает вопрос: в каком смысле расходящийся ряд (1) 
представляст функцию (3)? Чтобы ответить на этот вопрос, 
заметим, что при О, 1, 2) ws 
т. 


i. м (— ху +1 
Гм ey {— #9" -- T+ xt 


n=0 


и, следовательно, 


И = 5» (<) + К» <, (4) 
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где 
= МЛ" ° ый ° (5) 


п=0 


— частичная сумма ряда (1), а 


no 
етич 


Ки (X) = (xe {or (6) 


0 


— ero остаточный член. 
Если Вех = 0, то 


Г-Н ХЕ [71 = 1 


Аа т-- Ох", rae 


С другой стороны, если Rex < 0, у=агы хит < ут, 
то 
ПН хЁ|”' = | созее|, 


[Аи (20) | << (m1) tf 


В обоих случаях остаточный член имеет тот же порядок, 
что и первый отброшенный член ряда S (x), и быстро стре- 
мится к нулю, когда х—*0. Сходимость равномерна в любом 
секторе | ага х| == п— в, => 0. Если Re x = 0, то остаточный 
член меньше по модулю первого отброшенного члена, причем 
если x >>0, то остаток имеет тот же знак, что и первый 
отброшенный член. Таким образом, при x >>O0 ряд (1) ведет 
себя подобно схолящимся знакочередующимся рядам, за исклю- 
чением того, что остаточный член не стремится к нулю. 
Минимальное значение погрешности определяется наименьшим 
членом ряда (1) (получающимся, когда m примерно равно x). 

Теория асимптотических рядов ведет начало от Стил- 
тьеса (1886) и Пуанкаре (1886). Эту тсорию можно разбить 
на две части. В одной из них, теорни асимптотических рядов, 
изучаются такие вопросы, как «суммы» асимптотических 
рядов (‹«асимитотические пределы», «асимитотическая сходи- 
мость») и операции над асимитотическими рядами (алгебраи- 
ческие операции, дифференцирование, интегрирование, под- 


Pes cosec 


‚ где Вех< 0. (8) 
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становка асимптотических разложений переменного в асимпто- 
тические и сходящиеся ряды, содержащие это переменное, 
и тому подобное). Наиболее полное изложение этой теории 
можно найти в лекциях ван дер Корпута (1951, 1952) и его 
текущих работах. В этой книге мы ограничимся кратким вве- 
дением в теорию асимптотических рядов, обратив основное 
внимание на другую часть изучаемого вопроса — теорию 
асимптотических разложений. Центральной темой является 
здесь построение и изучение рядов, асимптотически пред- 
ставляющих данные функции. Эти функции часто задаются 
интегральными представлениями или степенными рядами либо 
возникают как решения дифференциальных уравнений; в пос- 
леднем случае переменное, по которому производится асимпто- 
тическое разложение, является либо независимым переменным 
в дифференциальном уравнении, либо входящим в это уравне- 
ние параметром. 
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ГЛАВА 1 
АСИМПТОТИЧЕСКИЕ РЯДЫ 


1.1. О-символика 


Как правило, под независимым переменным мы будем 
понимать вещественное или комплексное переменное; но в этой 
главе, если не оговорено противное, х означает переменный 
элемент топологического Т,-пространства (хаусдорфова про- 
странства). Обозначим через Ю множество значений, прини- 
маемых переменным х, а через x, некоторую предельную 
точку множества А (которая может и не принадлежать №). 
Символы 9(х), O(x) и аналогичные им обозначают функции 
от х, определенные на множестве А и принимающие веще- 
ственные или комплексные значения. 

Мы будем применять следующие обозначения. Если 
существует постоянная A (то есть число, не зависящее OT x), 
такая что || А|Ф| для всех x€R, мы будем писать 
о ==0($). Если существует такая постоянная А и окрестность 
U точки x,, что jy{<A|O] для всех x, принадлежащих 
ИПК, мы будем писать 9==O(%) при х—+х,. Аналогично, 
мы будем писать ф==0(ф) при x—x,, если для любого 
=> 0 найдется такая окрестность U, точки x,, что [9 |< =|%| 
для всех x, принадлежащих (И.П А. Если на множестве Ю 
функция ф не обращается в нуль, то эти три условия можно 
сформулировать проще: 9 = О ($) в В (при х—х, в А}, если 
ф/ф ограничено BR (при х—+х, в А); $=0(1) при x—x,, 
если g/b—+0 при x—+x, в А. Соотношения y= O ($) р 
ф==0(Ф) мы будем называть соотношениями порядка. 

В следующих примерах х является комплексным перемен- 
ным, а Sy — областью; 0 | х| < оо, | га х| < п/2 —А. 
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()e7* = O(x*), e7* == 0 (x"), если x—co B Ss, А>0, 
&« произвольно; эти соотношения не выполняются (при любом 
a), если A<O. 

(il)e"*==O(x’), если x—+0o в S,, при условии, что 
Кея 0; в области Rea < 0 это соотношение не имеет места. 

(Ш) е`*=0(х“) B Sy, если либо А > 0 и Вед <0, либо 
A=0 и Rea=0O. 

Читатель без труда докажет сформулированные выше 
утверждения. 

Всли функции, входящие в соотношения порядка, зависят 
от параметров, то, вообще говоря, постоянная A и окрестности U, 
U,, входящие в определения, зависят от этих параметров. 
Если A, U, О. могут быть выбраны независимо от парамет- 
ров, то говорят, что соотношения порядка выполняются 
равномерно относительно этих параметров. 

Операции ¢ соотношениями порядка выполняются по неко- 
торым простым правилам. Мы приведем наиболее часто встре- 
чающиеся правила для О-символа; соответствующие правила 
справедливы и для о-символа. В указанных ниже правилах R 
и x, фиксированы, и слова «при х—+х,» опускаются. 

Если ¢==O(%) их>0, то 


[97 =0( |4”. (1) 


Если 9; ==О(Ф;), {=1, 2, ..., А, и я; — постоянные, то 
Уи = О). (2) 
1 a 


Эти соотношения выполняются также для бесконечных рядов 
при условии, что $; = О(ф;) равномерно по 2. В случае бес- 
конечных рядов равенство (2) и аналогичные соотношения 
имеют следующий смысл. Если ряд У|9/5,| сходится, то 
сходится и ряд Уа, причем выполнено соотношение (2); 
если же ряд У [ай,| расходится, то о ряде Ук; ничего 
не утверждается. 

Если 9;==O(%;), i=1, 2, ..., А; а, постоянные и 
|4; |= при 1—1, 2, ..., А для всех х, принадлежащих 
пересечению множества А и некоторой окрестности И, 
точки xX,, то 


= O(Y). (3) 
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Это соотношение справедливо для бесконечных рядов прн ус- 
о 
ловии, что 9;==0О(ф‚) равномерно no i u У [а;| < оо. 
Если ¥,=O(%,), i=1, ..., А, то 


Пи=о И. (4) 


г 


Доказательство соотношения (1) очевидно. Для того 
чтобы доказать соотношение (2), возьмем числа А; и окрест- 
ности (; точки x,, соответствующие функциям ®,. Если число 
этих функций конечно, то существует число A, большее, 
чем все числа A; и окрестность U, содержащаяся во всех 
окрестностях U;. Тогда, если x принадлежит (ИП R, выпол- 
няется неравенство 


[Узи |= У 


чем и доказано соотношение (2). Если множество функций 9, 
бесконечно, то существование A и U вытекает из равномер- 
ности соотношений порядка по i. Соотношение (3) может 
быть выведено из (2), так как при сделанных предположениях 
можно выбрать окрестность U так, чтобы она принадлежала 
и U,, a тогда 


АИФ АУ а: ||4, 


АУ: ||: = А [= А, 


где А, =АУ|а,| — постоянная. Доказательство (4) анало- 
гично доказательству (2). 

Соотношения порядка можно интегрировать как по неза- 
висимому переменному, так и по параметрам. Ради простоты 
ограничимся интегралами по вещественным переменным, однако 
полученные результаты можно распространить на комплекс- 
ные и абстрактные переменные. 

Пусть х-— вещественное — переменное, —Ю — интервал 
а<х<Ь, из =0(ф) при х —6. Если функции ди ф из- 
меримы в R, то 


b b 


J oat=O(§ [40 при х—6. (5) 
* x 
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ь 


Доказательство. Если ре [а = оо, то утверж- 


х 
b 


дение очевидно. Если ра < во для некоторого x, то 


& 6 
существуют такие числа А и X, что 
b ‘ 

оо и |9) = 44| apn Xx <p, 
x 
а тогда 
b b b 
а = = (|4 при Х<х<Ь. 
x x x 
Пусть х — переменный элемент множества R в хаусдор- 
фовом пространстве, у — вещественный параметр, a<y <8 
иф(х, y)==O(4(x, у)) равномерно по у, когда х —х,. Если 
для любого фиксированного x в Ю gud являются измеримыми 
функциями от у в интервале о < у<В, то 
8 4 
Jets yay Ol (а, УЧ при xx (6) 
a a 
Доказательство аналогично доказательству соотношения 
(5). Пользуясь равномерностью О-символа, можно выбрать А 
и Ц, не зависящие от у так, что || < А|%|. Интегрируя 
это неравенство по у, мы приходим к соотношению (6). 
Вообще говоря, недопустимо дифференцирование соотно- 
шений порядка как по независимому переменному, так и по 
параметрам. Однако некоторые общие результаты о диффе- 
ренцировании соотношений порядка имеют место в случае 
аналитических функций комплексного переменного (см. п. 1.6). 
В заключение этого пункта приведем некоторые формулы, 
касающиеся комбинирования соотношений порядка: 


0(0(¢)) =0(9), (7) 
0 (0 (9) 0 (0 (9) = 0 (0 (d)) = 0 (9), (8) 
O(~yO(h) = Oleh), ` (9) 
Olgyo(d) =0(4)0(b) =o ($9), (10) 
О-о = O(4) +0 (4) = O(%), (11) 


0(¥)-+-0(~) ==0(4). (12) 
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Эти очевидные формулы могут быть распространены на ком- 
бинации любого конечного числа О- и о-символов. 

При ссылках на приведенные правила мы будем указывать 
номер соотношения, выражающего окончательный вывод, 
причем тот же самый номер будет употребляться для соот- 
ветствующего правила относительно о0-символа. Например, 
соотношение (1) указывает либо правило, что из © ==О ($) и 
а > 0 вытекает | $ "== 0 (|ф |”), либо правило, что из 9 ==0 ($) 
uw a>>0 вытекает |9’ ==0(|$|”). 


1.2. Асимптотические последовательности 


В этом пункте А, x, X,, ф имеют тот же самый смысл, 
что и в п. 1.1. Конечная или бесконечная последователь- 
ность функций 9, 9,, ...’ сокращенно будет обозначать- 
ся {9}. 

Последовательность функций {o,} называется асимпто- 
тической последовательностью при х-—>х, в К, если для 
всех п функции 4, определены в Ки ¥,,,=0(%,) при 
хх, в В. 

‚Если последовательность  бесконечна и $,.,=0(9,) 
равномерно но 2, то {2„! называется асимптотической после- 
довательностью, равномерной по п. Если ф„ зависит от пара- 
метров и 9, ,,==0(%,) равномерно относительно параметров, 
то {9„} называется асимптотической последовательностью, 
равномерной по параметрам. 

Приведем некоторые примеры асимитотических последова- 
тельностей. В этих примерах х — комплексное переменное; 
R, за исключением особо оговоренных случаев, — комплексная 
плоскость и 5. — область, определенная в п. 1.1. 


(о {(«—-x,)"}, хх 
(II) 4х" x—x,; 
(Ш) {Хх}, хм B Sy, 
где Вел, ., > Ве), для всех 7; 
(У) {х-*} хо, 
где h, вещественно и },,, >», для любого п; 
У ея, 
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где хи», либо обладают теми же свойствами, что ив (III), 
либо теми же свойствами, что в (IV); 


(VI) {2-75 


где x—co в Sy, и либо А220, причем Ве%,,, ->Rel, 
для всех и, либо А > 0 и )., произвольны; 


(VL) Гы Pw-+a}, хм в Ss, A>—n/2 


Читателю полезно проверить, что последовательности 
(1) — (УП) являются асимптотическими, и выяснить смысл огра- 
ничений, наложенных в этих примерах на А и *,. Например, 
почему последовательность (Ш) при произвольном п), не 
является асимптотической, когда х — со в комплексной пло- 
скости (но является ре oat в 5.)? Бесконечная 
последовательность {Г (х— и) Г(х--п)}, п=1, 2, ..., не 
является асимптотической при X— со в любой области, содер- 
жащей неограниченную часть вещественной. оси; однако она 
является асимптотической при х —+ со в любой области, замы- 
кание которой целиком лежит в верхней или в нижней полу- 
плоскости. Конечная последовательность {Г (х — п)/Г (х-Ёп)}, 
n==1,2,..., М, является асимитотической при х —> oO в лю- 
бой области К. 

Если заданы некоторые асимптотические последователь- 
ности, то можно получить новые асимптотические последова- 
тельности с помощью процессов, опирающихся в основном на 
приведенные в п. 1.1 правила действий над О- и о-симво- 
лами. При описании некоторых из этих процессов мы ограничимся 
вещественными переменными, однако их можно распространить 
и на переменные более общего характера. В большинстве 
случаев x, и R не упоминаются; в этих случаях OHH пред- 
полагаются фиксированными. 

„Любая подпоследовательность асимптотической после- 
довательности является асимптотической последователь- 
ностью. Доказательство вытекает из 1.1 (8). 

Если 4%9,} — асимптотическая последовательность и 
a >>0, mo {| ,[} также является асимптотической после- 
довательностью. Доказательство вытекает из 1.1(1). 

Две последовательности {9„} и {b,}, такие, что ф,==0О($,) 
и Оф») при всех 7, называются эквивалентными. Если 
{y,} и 4,} — эквивалентные последовательности, причем 
р Их последовательность, то \%„} так- 
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же является асимптотической последовательностью. Для 
доказательства заметим, что, в силу 1.1(8), имеем 


Фин: = О (9,1) = O (9 (9,)} = O(0 (0 (ф,))) = 0 (4,)- 


Если {9} и {,} — асимптотические последовательно- 
сти, имеющие одно и то же число членов, mo 44,b,} 
является асимптотической последовательностью. Доказа- 
тельство вытекает из 1.1(10). 

Пусть {9}, П=,..., М, — асимптотическая последо- 
вательность ua, „П==1,..., М; 1=0,1,..., №, # < М,— 
множество положительных чисел, таких, ито 9... = а 
при всех п, i; если 


п 


k 
= til Sarl Hl, NR, (1) 


(=o 


mo {,\ является асимптотической последовательностью. 
В этом утверждении N может быть конечным или бесконеч- 
ным, в то время как А конечно. Для доказательства утверж- 
дения заметим, что, в силу конечности А для любых пи => 0, 
найдется окрестность (/, точки x,, такая, что при г==л, п- 1, 
wee, A+R в О.П Ю выполняется неравенство |y,,,|< 
<<(¢,|- aa 


|=, (2) 


lou: 
Oni Ynti 


Gna =. нии Pnsini | 5 
ПЕ 


Следующая теорема дает распространение этого утверж- 
дения на бесконечные ряды. 

Пусть \9„} равномерно по п является асимптотиче- 
ской последовательностью и пусть ор И i ee 
1=0,1,..., — множество положительных чисел, таких, 
что O41; S%,; для всех пи 1. Положим 


ea 

„= оно МНЕ aes (3) 
Если бесконечный ряд, определяющий %,, сходится в неко- 
торой окрестности точки X,, то существует такое под- 
множество R, в R, для которого x, является предель- 
ной точкой, что ряд (3) сходится в R, и {,) равномер- 
но по п является асимптотической последовательностью 
при х-—х, в R,. 
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Доказательство. Так как последовательность {9,} 
является асимптотической равномерно по п, то существует 
подмножество А, в К, такое, что x, является предельной 
точкой для К, и |=, | для всех x из R, и всех п. 
Если x принадлежит R,, то 


Dens tl Фин = vad ete | < У, |4 Sisal 


и потому бесконечные ряды (3) мажорируются рядом для Ф,. 
По условию ряд для \, сходится в некоторой окрестности R, 
точки x,. Положим А, =А, П R,. Все функции Ф, опреде- 
лены в Ю,, х, является предельной точкой для R, и, поскольку 
{9,} являстся асимптотической последовательностью равно- 
мерно по м, равенство (2) при Ё—=со выполняется равно- 
мерно по 7. к 

Новые асимитотические последовательности могут быть 
получены © помощью интегрирования двумя различными пу- 
тями. 

Если {$„(х, y)} равномерно по у, «< у<В, является 
асимптотической последовательностью при х—х, в Ки 
если сходятся интегралы 


3 
®, (x)= |, ey) lay, 4) 


то {®,} является асимптотической последовательностью. 
Доказательство вытекает из 1.1(6). Как и в случае (3), 
достаточно допустить, что $, (х, у) являются измеримыми функ- 
циями по у и что 9, — суммируемая функция. Тогда из того, 
что интеграл {lent ay мажорируется интегралом Vie, tay, 
вытекает, что все функции Y, также суммируемы, быть может, 
на некотором меньшем множестве R,. 

Если х-— вещественное переменное, КЮ — интервал 
а<х<Ь, {1} — асимптотическая последовательность при 


х—Ь и если интегралы 
b 


®, (x) =, (flat (5) 


сходятся, то {Ф„} является асимптотической последова- 
тельностью при x —> b.- Доказательство вытекает из 1.1 (5), 
причем, как и выше, достаточно предположить, что все функ- 
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ции ©, измеримы, а функция 9, суммируема; в этом случае 
результат имеет место на некотором интервале а, < х<5. 

Заметим, что дифференцирование асимптотических после- 
довательностей не всегда приводит к асимптотической после- 
довательности. Например, положим 


$,= x "[a-+cos(x”)[, п==1,2,... 


Тогда {y,} является асимптотической последовательностью, 
когда х—› со на вещественной оси, но {g’,} не является асимп- 
тотической последовательностью. 


1.3. Асимптотические разложения 


В этом и следующих пунктах x, х., R имеют тот же самый 
смысл, что иви. 1.1. Через {9,}, {dat {yn}... всегда 
обозначаются асимптотические последовательности при X —> x, 
в А; f(x), g(x), A(x), ... — функции от x, определенные в Ю 
и принимающие числовые значения; а, b, с,... — постоян- 
ные (то есть не зависят от X). 

Формальный ряд У, (x) называется асимптотиче- 
ским разложением функции f(x) при х—х, Co М№-го члена 
(включительно), если 

м 


F(x)= У а, (x) olp,) при x x,. (1) 


n=) 
Мы будем обычно обозначать асимптотическое разложе- 
ние до №-го члена следующим образом: 

f(x) — Уаз, (®) до М№-го члена, когда x—+x, BR, (2) 
опуская, как правило, слова «в R». Асимптотическое разло- 
жение, состоящее из одного члена, записывается в виде 

Л(х) > ag, (x), х-х, (3) 


и называется асимптотическим представлением. Асимптоти- 
ческое разложение, справедливое до любого члена (то есть 
при №М==со), записывается в виде 


19 > Хан, (x), хх (4) 


и называется просто асиймптотическим разложением. Асимп- 
тотическое разложение может как сходиться, так и расхо- 
диться. В большинстве книг рассматриваются только случаи 
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N=1 и Noo, но мы будем считать N любым целым поло- 
жительным числом. ’ 

Если асимптотическое разложение до N-ro члена, где № 
конечно, зависит от некоторых параметров, мы будем говорить, 
что оно выполняется равномерно по этим параметрам, если 
остаточный член в формуле (1) является o (7х) равномерно 
по параметрам. Асимитотическое разложение (№ == со), зави- 


сящее от некоторых параметров, называется выполняющимся 
м 


равномерно по этим параметрам, если f— > 0, Vy =0(9м) 
Re 

равномерно по параметрам для всех достаточно больших зна- 

чений М (но не обязательно равномерно по M). 

Формальные (конечные или бесконечные) ряды Уаз, 
называются асимптотическими рядами. Вели ~,—=x~'", 
то говорят OO асимптотическом ряде степеней, а если 
9„=<х+”, то 06 асимитотическом стеленндлм ряде. Например, 
Уши! (—х)"`' является  асимитотическим  разложе- 
нием в степенной ряд ups х—*0 в области Son +: функ- 
ции f(x), определенной равенством (3) во ввелении. Некото- 
рые авторы говорят об асимитотических стеиенных рядах, если 
Pp ==, (x) xt", но в этом случае удобнее называгь асимпто- 
тическим степенным рядом такой ряд, деленный на 9, (x). 

Из равенства (1) вытекает, что коэффициенты асимптоти- 
ческого разложения до N-ro члена могут быть вычислены 
с помощью следующей рекуррентной формулы 


lim {И (x) — юз 29 (X)]]9q (xp, т=1,..., М. (5) 


X> ль n=1 


Обратно *), предположим, что мы имеем №--1 функций, 


1х), 9, (х), 00s Gy (%)s 


определенных в А. Если выполнено равенство (5) и @,„=20 
при т==1,..., №, mo {9,} является асимптотической 
последовательностью при х—х, и > a,9, — асимптотиче- 
ское разложение Oo №-го члена функции f(x) при х-ъх,. 

Чтобы доказать, что 12} является асимитотической по- 
следовательностью, надо показать, что %,,4,==0(%,) при 


`_ *) Эта теорема принадлежит Макки (A. G. Mackie). 
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m=, Pts ‚ N—1. Ho из (5) вытекает 


Сравнивая последние два равенства, получаем соотношение 
[1-70 (1) Зин = 9 (9). 


Если @,4,520, то a,,,-0(1) отлично от нуля для всех x 
В некоторой окрестности xy. Мы можем поэтому разделить на 
это выражение и получить 9„., ==0(2„). Таким образом, {2 „} 
является асимптотической последовательностью. Кроме того, 
полагая в равенстве (5) т==М, убеждаемся в выполнении 
соотношения (1), а потому У, является асимитотическим 
разложением до N-ro члена функции /. 
Если Уаз, (х) является асимитотическим разложением до 
№-го члена функции f(x), то тот же самый формальный ряд 
является асимитотическим р тожением этой функции и до 
М-го члена ММ. Мы имеем, таким образом, несколько 
более точный результат 
м 

f(x) == Dag, ©) + Oley as) хх, M=1,...,N—1, (6) 
n=1 

являющийся непосредственным следствием соотношения (1). 

Если x, и Ю фиксированы, то из равенства (5) вытекает, 
что при заданной асимитотической последовательности асимп- 
тотическое разложение данной функции, имеющее заданное 
число членов, однозначно определено. С другой стороны, одна 
и та же функция может иметь асимптотические разложения 
по двум различным ACHMHTOTHYeCKHM последовательностям, при- 
чем эти последовательности могут не быть эквивалентными 
в смысле п. 1.2. Например, 


| м f\tetyon № 
тя у (— 1)" 'x7", x —+ 00, 

1 x an 
Ears (x—1)x7", хо, 


1 co ae зп 
ret DIO et pt, x00, 
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В этом примере все три асимптотических разложения являются 
сходящимися рядами, когда |x|>>1. Часто бывает, что неко- 
торые асимптотические разложения функции расходятся, в то 
время как другие сходятся. Преобразование расходящихся 
асимптотических разложений в сходящиеся интересно с анали- 
тической точки зрения, но мало дает для вычислительной 
практики. Преобразование асимптотических разложений в схо- 
дящиеся разложения, а также в разложения, более удобные 
для численных расчетов, изучались Эйри (1937), ван дер 
Kopuytom (1951), Миллером (1952), ван Вингаарденом (1953), 
Ватсоном (19125) и др. 

Асимптотическое разложение не определяет однозначно 
суммы f(x). Например, функции (1 x)7", (1 +e7*)/(1 +x), 
(ite Ух 1-х)-' имеют одно и то же асимптотическое раз- 
ложение У (—1)"'х7", когда х-+ оо в Ss, A>0. Если 
задана (конечная или бесконечная) асимптотическая последо- 
вательность {9,} при х—+х, в К, то множество всех функ- 
ций, определенных в Л, распадается на классы попарно экви- 
валентных функций. Функции f(x) и g(x) называются асимл- 
тотически равными относительно {2}, если для всех п 
выполняется соотношение 


F(x) — &(х) ==0(%,) при x—>x, в К. 


Поэтому асимптотический ряд представляет не одну функцию, 
а класс асимптотически равных функций. 


1.4. Линейные операции над асимптотическими 
разложениями 


Если F~ Ув 00 М№-го члена и g~ УВ, д0 М-го 
члена и если &, В — некоторые числа, то 


af (x)-+ Bg (x) ~ У (2a, 86) ф,(х) д0 N-20 члена. (1) 


Эта теорема очевидна, равно как и ее обобщение на линей- 
ные комбинации любого конечного числа асимптотических раз- 
ложений. Обобщение этой теоремы на бесконечные ряды, 
состоящие из асимптотических разложений, формулируется 
следующим образом: 

Пусть равномерно по i, 1=1,2,..., имеем f;(x)~ 
~ Deni Fn (x) до №-го члена; пусть, далее, числа a; таковы, 
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a 
что ряд > а; абсолютно сходится, а ряды 


ice} 
An = 2 а, Ва (2) 
i= 
сходятся при всех п. Тогда ряд das (x) cxodumca 6 He- 
которой окрестности точки x,, причем 
[2 


Р(х) = > aif) (x) — У Ане, (x) д0 N-20 члена. (3) 


f=? 


Доказательство. Мы имеем, равномерно no в 


A= 


fir, ‚бы, n= 0 (¥y)- 


i 


Кроме того, })|2;|<co. В силу 1.1(3), 
оо N 
>>} 2 (Fi х 2, 4.) =0 ($), 
ti=1 пы! 
причем стоящий слева ряд сходится в некоторой окрестности 
точки x,. Всли N< oo, то, прибавляя к обеим частям этого 
№ 


равенства У A,%,, мы получаем выражение (3). Если № == со, 
n=1 

то как предположения, Tak и заключение теоремы справед- 

ливо для всех достаточно больших №, а потому У А, 

является асимптотическим разложением функция F(x). 

Формулу (1) можно обобщить на конечные или бесконеч- 
ные асимптотические разложения. 

Пусть {9}, п=1, ..., N<oo, и {b,}, т=1,... 
sees М = оо, являются асимптотическими послебователь- 
ностями для обних и тех же Ru x,, причем Oy=O(d,) 
при любом т; если f~ Хану, до N-2o члена и, при любом п, 
4, ^ УВ,» д0 М-го члена, то 


F(X) > Ус,» (x) до М-го члена, (4) 
где 
м 
C= У аби. (5) 
a 
Пусть {9,}, 2n=1,2,..., и {d,}, m=1,..., M<oo, 


являются асимптотическими последовательностями, при- 
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цем для л0бого п найдется такое целое число в (п) = М, 
что РА, gotta пс и 9, =0(, (и); если 
{~ da, Qn» On У, nim 90 М-го члена равномерно по п, 
ряд ый абсолютно ‘сходится и бесконечные ряды (5) 
сходятся при любом т, то имеет место равенство (4). 

Доказательство очевидно, если M и М — конечные числа, 
так как тогда, в силу 1.1(3), 


уу == 3 ао (Sy) = 
N 


1 
M 


=D aD bide Роб + 0 (at) = 


т=! 


м 
= У сом). 


ma) 


Если M==00, то это же самое рассуждение справедливо 
при любом М и, следовательно, deb, является асимито- 
тическим разложением функции f(x). Для обобщения на слу- 
чай № = со надо использовать обобщение формулы 1.1(3) 
на бесконечные ряды. 

Перейдем теперь к вопросу об интегрировании асимито- 
тических разложений как по вещественному параметру у, 
так и по переменному х; в последнем случае мы предпола- 
гаем, что переменное х вещественно, 

Пусть f(x, у) > a; (уф, (x) до N-20 члена равномерно 
по у, a<y<cB, пусть, далее, f(x, у) при любом фикси- 
рованном x и а„(у), при любом п, являются измеримыми 
функциями от у, а hly) является суммируемой функцией, 
такой, что все интегралы 


А, = J hy), (y) dy (6) 


сходятся. Тогда Can всех х в некоторой окрестности 
точки X, сходятся интегралы 
3 
Ре) = в (У), y) dy, (7) 
причем 
F(x) ~ У Ар, (х) 90 М-го члена. | (8) 
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Эта теорема доказывается почти так же, как и соответ- 
ствующая теорема для бесконечных рядов, за исключением 
того, что вместо 1.1(3) надо использовать 1.1(6). Очевидно, 
что в этой теореме интервал (х, в можно заменить любым 
измеримым множеством конечной или бесконечной меры или 
множеством в многомерном пространстве. 

Пусть х— вещественное переменное, КЮ — интервал 
ALKA, а {3 (x)} — асимптотическая последовательность 
положительных функций при х—65, такая, что все интес- 
ралы . 


(9) 


Ф,„ (= 


ес 
5 
x 
as 
ae 
& 


схобятся. Если Г(х)— Уа„т,(х) 00 М№-го члена при х — 5 
и /(х) — измеримая функция, то интеграл 


b 
F(iy=\ shat (10) 


сходится для всех х, лежащих в некотором интервале 
с3х<Ь, прием 


F(x) > Уа,Ф, (x) Co М-го члена, когда x—>b. (11). 


Доказательство вытекает из 1.1(5). 

Лифференцирование асимптотических разложений как по 
переменному xX, так и по параметру, вообще говоря, нело- 
пустимо. Некоторые общие результаты относительно дифе е- 
ренцирования асимитотических `разложений аналитических 
функций комплексного переменного будут указаны в п. 1.6. 


1.5. Другие огерации над асимптотическими 
разложениями 


Умножение асимптотических рядов не приводит, вообще 
говоря, к асимптотическому ряду, так как в формальное 
произведение рядов Уа,р, и УЬ,2, входят все произведе- 
ния 9,9, а вообще говоря, невозможно упорядочить систему 
функций {9,0}, т, п=1,..., М, так, чтобы получилась 
асимптотическая последовательность. В некоторых важных 
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случаях, однако, либо произведения 9,0, образуют асимпто- 
тическую последовательность, либо эти произведения обла- 
дают асимптотическими разложениями по некоторой асимпто- 
тической последовательности (не обязательно совпадающей 
с {%,}). Докажем сначала общий результат об умножении 
двух асимитотических разложений. 


Пусть {gph == Nye. М, (Yah MEL, ves М, и {д}, 
k==1,..., А, являются такими, асимитотическими последова- 
тельностями, что м == О(ук), 9, = О (ук) и 

Фи — Di Cnmete ДО К-го члена. (1) 


Если f~ У ан, _д0 N-20 члена и У, G0 М-го 
члена, то fg~ УС,], 00 K-20 члена, где 
м 


м 
C= 2 DPC nme: (2) 


n= ms 


Здесь К может быть конечным или бесконечным, aNuM 
конечны. Результат остается справедливым, если М или М, 
или Ми М одновременно бесконечны при условии, что все 
бесконечные ряды (или двойные ряды) в формуле (2) схо- 
дятся. 

Коэффициенты С, получаются при умножении рядов Ув 
и У,» я подстановке разложения (1), поэтому вместо 
равенства (2) можно сказать словами: «где коэффициенты С, 
получаются с помощью формальной подстановки». Это выра- 
жение будет употребляться далее в аналогичных случаях. 

Докажем сначала теорему для конечных №, М, К. В силу 
соотношений (1) и (2), имеем 


у N м 
ПХ ая, Но ИХ в, to Call = 


м М 
= У У аб --о (2) + ig (gy) = 


к mot. 
= = Cee 7 9 (Lx) офи) + obs)» 


чем требуемый результат и доказан. Если разложение (1) 
справедливо до любого члена, причем о Фи и Yad, являются 
O(y,) для любого k, то проведенное выше вычисление 
сохраняет силу для любого К, а потому полученный резуль- 
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тат справедлив и при К== оо. Обобщение на случай беско- 
нечных Ми М доказывается подобным образом при условии, 
что ряды, определяющие коэффициенты C,, сходятся. 

Последовательность функций {9„}, п=1,..., М, назы- 
вается мультипликативной асимптотической последова- 
тельностью, если {,} является асимптотической после- 
довательностью, 9, =0(1) и 9. ~ DY Cume2e 90 N-20 
члена, т, п=1,2,..., М. В случае мультипликативных 
асимптотических последовательностей полученный выше ре- 
зультат об умножении асимитотических разложений может 
быть существенно усилен. 

Если {9,}, П==1,..., М, является мультипликативной 
асимптотической последовательностью, 


Л - Zan, 4, д0 N-20 члена, i=1,..., №, 


u P(2,, ..., Zp) является многочленом от Е комплексных 
переменных 2,..., 2» то функция F(x\==P(f,, ..., fy) 
имеет асимптотическое разложение Е — У А, G0 №-го 
члена, коэффициенты которого А, могут быть вычислены 
с помощью формальной подстановки. 

Для того чтобы доказать эту теорему, заметим, что 
в случае мультипликативной асимптотической последователь- 
ности {9„} мы имеем 2 9,==0О(9л), а функции 9,9, имеют 
асимптотические разложения до №-го члена. По общей тео- 
реме отсюда вытекает, что если 1- Уаз до N-ro члена 
и > Ув, до N-ro члена, то функция fg обладает асимп- 
тотическим разложением Donen до N-ro члена, коэффици- 
енты которого могут быть вычислены с помощью формальной 
подстановки. Вычисление любого многочлена P(f,, ..., fy) 
сводится к конечному числу операций, каждая из которых 
заключается либо в образовании линейной комбинации, либо 
в умножении двух асимптотических последовательностей. 
Каждая из этих операций сохраняет асимптотический харак- 
тер разложения, причем получающиеся при этом разложения 
могут быть вычислены с помощью формальной подстановки. 
Тем самым наша теорема доказана. 

Полученный для многочленов результат может быть при 
некоторых предположениях распространен на (сходящиеся} 
степенные ряды и даже на асимптотические степенные ряды. 
Ради простоты ограничимся случаем одного переменного 2, 
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хотя соответствующее обобщение этого результата справед- 
ливо и в случае многих переменных. 


Пусть {2,}, п=1,..., М, является мультипликативной 
асимитотической последовательностью, такой, что $, =0(1) 
и м = O(9y) для некоторого натурального М. Если 


Г(2) > Ус” 00 М-го члена, когда z—+O в комплекс- 
ной плоскостя и 2=2(х)^ Day, 00 М-г0 члена, 
когда х-—>х, в R, то функция Р(х) = (2 (х))} обладает 
асииптотическим раз. 1ожением Уд, do N-2o члена, 
когда X—>X,, причем коэффициенты A, этого разложе- 
ния могут быть вычислены с помощью формальной под- 
становки. 
Доказательство. Из сделанных предположений выте- 
. KaeT, что функция ze имеют асимптотические разложения 
=" — VO gn G_, до №-го члена, причем z@==O(9,). Следова- 
тельно, можно применить теорему из п. 1.4 о подстановке 
одного асимптотического разложения в другое асимптотическое 
разложение. Теорема доказана. 
Рассмотрим теперь функции f(x), имеющие асимптоти- 
ческое разложение вида , 


N 

Пе) e+ Хар, +0 ($y). (3) 
В этом случае из теоремы вытекает, что если с = 0 и функ- 
ции {9,} удовлетворяют условиям теоремы, то [/(x)]~* 
обладает асимитотическим разложением того же вида. Иными 
словами, асимптотические разложения вида (3) можно делить 
друг на друга. Это позволяет распространить полученные 
выше результаты на рациональные функции, 

Ecau {9 П=Т,..., М, является мультипликативной 
асимптотической последовательностью, 2, =0(1), |, |= 
=0(7м) для некоторого М, f;~ Den, Pn 00 N-20 члена, 
1=1 ..., № и Plz, ..., 2) является рациональной 
функцией от Е комплексных переменных 2,,..., Zp, зна- 
‚менатель которой отличен от нуля при 2.==2,==... 
oe ==2,==0, mo функция Е(х)==Р(,...,/,) обладает 
асимптотическим разложением A,-+- У A,2, до №-го члена, 
коэффициенты A, которого могут быть вычислены с по- 
мощью формальной подстановки. 

При тех же самых предположениях существует асимпто- 
тическое разложение функции g(F(x)), где g(t) является 
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функцией комплексного переменного €, регулярной в некото- 
рой окрестности точки &, = Р (0, ..., 0). Таким путем можно 
доказать существование асимптотических разложений для 
таких выражений, как, например, ехр [Р(Х,, ..., /»]|. 


1.6. Асимптотические степенные ряды 


Последовательность функций {х-"}, п=0, 1, 2, ..., 
или n==1, 9,..., является мультипликативной асимитотиче- 
ской последовательностью при х—> сою в любой области 
комплексной плоскости, не содержащей точки х==0. Эта 
последовательность удовлетворяет всем условиям, наложенным 
на асимптотические последовательности в двух предыдущих 
пунктах, за исключением того, что в некоторых теоремах 
п. 1.5 следует исключить значение п==0. Кроме того, она 
обладает некоторыми специальными свойствами. 

Асимптотическое разложение 


f(x) ~ a, ++ ие -- ие ... до N-ro члена, когда x — со, (1) 


мы называем асимитотическим степенным рядом. Из pe- 
зультатов пп. 1.5 и 1.6 вытекает, что асимптотическое 
разложение в степенной ряд можно умножать на постоянную 
величину, а два таких разложения можно складывать и умно- 
жать друг на друга, а также делить при условии, что в раз- 
ложении знаменателя @,>5£0. Асимптотические степенные 
ряды можно подставлять в конечные линейные комбинации, 
многочлены и рациональные функции при условии, что зна- 
менатель этих функций не стремится к нулю при х — oo; 
их можно подставлять также в сходящиеся и асимптотические 
степенные ряды D2" z—>0, при условии, что в разложе- 
нии (1) функции z== f(x) коэффициент а, =0. Подстановка 
разложения (1) в асимптотические и сходящиеся ряды другого 
типа возможна, если выполнены условия, указанные в п. 1.5. 
Во всех этих случаях коэффициенты новых разложений полу- 
чаются с помощью формальных подстановок и перестановки 
членов. Если асимптотическое разложение в степенной ряд (1) 
имеет место равномерно по некоторому параметру, то его 
можно интегрировать по этому параметру. Наконец, если 
имеет место разложение (1), то функция Л(х) — а —а,|х 
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интегрируема, причем 


в [po—a—4] at~ 24 H+ eet... @) 


до (N— 2)-го члена, когда х— oo. 

Простым следствием последнего результата является сле- 
дующая теорема о дифференцировании таких разложений. 
Если функция f(x) в равенстве (1) дифференцируема, 
а функция Г’(х) может быть разложена в асимптотине- 
ский степенной ряд, то 


л®- 


00 N— 1-20 члена, когда x—+ со. (3) 


В случае аналитических функций имеет место более 
точное утверждение, именно, можно не предполагать, что 
функция /’(х) имеет разложение в асимитотический степенной 
ряд. Пусть К — область 

Ге, а«авх< в. 
Пусть, далее, a, >a, «я, < В, «Ви Ю, — область 
Jx|>>a,, а Зах В. 
Если функция f(x) регулярна в области R и разложение (1) 
имеет место равномерно по агох, когда x-+co в Ю, то 
разложение (3) имеет место равномерно по агох, когда 
х—+ о в К,. Доказательство этой теоремы вытекает из 
интегральной формулы Коши для производной функции f(x) 
1 1 (=) 
ы т 
ба. (4) 
Е 
При заданных К, А, найдется такое ¢>>0, что для любого 
x BR, окружность с центром в х и радиусом =|х! лежит 
в Ю; мы можем выбрать эту окружность в качестве контура 
интегрирования в формуле (4). Вдоль этой окружности 
z=x+exe* и 0 = {= 2m, а потому формула (4) принимает 
вид 


2n 


e~ Fix (1 +-ee")] dt. (5) 
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Но функция е`ИУ[(х-вхе")] имеет асимптотическое разло- 
жение в степенной' ряд равномерно по &, которое можно про- 
интегрировать по # отсюда вытекает, что функция f” (x) 
имеет асимитотическое разложение в степенной ряд. Как 
было показано выше, это разложение задается формулой (3). 

Аналитические функции могут, вообще говоря, иметь 
в разных секторах различные асимптотические разложения 
в степенные ряды (явление Стокса). Это не имеет места, 
если аналитическая функция регулярна на бесконечности, что 
вытекает из следующей теоремы. 

Если функция f(x) однозначна и регулярна при | х |`> а, 
причем формула (1) справедлива для всех значений arg x, 
то степенной ряд в формуле (1) сходится при достаточно 
больших значениях |x| и сумма его равна f(x). Чтобы 
доказать это, положим х==1[& и 8(0)=а,, 8(8) =/ (1/8), 
0< Е! <1|а]*. Тогда функция g(§) является однозначной 
непрерывной функцией в области |{|<|@|”', причем эта 
функция регулярна в этой области, за исключением, быть 
может, точки &==0. Но в точке Ё==0 функция g не имеет 
ни полюса, ни существенной особенности, так как она огра- 
ничена в некоторой окрестности этой точки, Таким образом, 
функция g(&) регулярна при &==0 и может быть разложена 
в ряд Маклорена. Из теоремы единственности для асимпто- 
тических разложений вытекает, что разложение (1) при 
х==1|$ является разложением Маклорена. 


1.7. Суммирование асимптотических рядов 


В п. 1.3 было отмечено, что асимптотическая последова- 
тельность {9„} задает соотношение эквивалентности между 
функциями, определенными в области Ю. Две функции, опре- 
деленные в этой области, асимптотически равны, если их 
разность для всех и есть 0(9,). Асимптотически равные 
функции имеют одинаковые асимптотические разложения. 
Поэтому, если задано асимптотическое разложение f — Ура», 
то естественно назвать суммой асимптотического ряда 
У Ф„ класс всех функций, асимптотически равных f. 

Мы закончим эту главу доказательством того, что любой 
асимптотический ряд имеет сумму. Результаты этого вида 
были доказаны для асимптотических степенных рядов Боре- 
лем и Карлеманом (1926), для рядов, мажорируемых асимпто- 
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тическим рядом степеней, — ван дер Корпутом (19546) и для 
асимптотических рядов аналитических функций — Карлеманом 
(1926). Излагаемое ниже доказательство является видоизме- 
нением доказательства ван дер Корпута. 

Асимптотический ряд является формальным конечным или 
бесконечным рядом Уа„9,(х), где {„} — асимптотическая 
последовательность и а, — постоянные. Так как подпоследо- 
вательность асимптотической последовательности является по- 
следовательностью того же вида, мы можем считать, что 
а,-20 при всех п. Асимптотическая сумма ряда Ув, яв- 
ляется классом асимптотически равных функций, и мы дока- 
жем существование асимптотической суммы, построив один 
из представителей этого класса. Если \@,Ф„ является ко- 
нечным асимптотическим рядом, то в качестве представителя 
асимптотической суммы можно взять сумму 


ag, + CL Нам 
понимаемую в обычном смысле. Поэтому достаточно дать до- 
казательство для бесконечных асимптотических рядов Ура, 
в которых а„5Е0 при всех п. 
Пусть U, — окрестность точки х, и пусть для всех л=1, 
2,... U, является такой окрестностью точки X,, что замы- 
Kanne U,,, лежит в О, и 


1 
[ан бич |= зари | 


для всех х, принадлежащих пересечению И, и Ю; такие ок- 
рестности существуют, так как 


Чин =O (Ayn) 
Обозначим для любого п через pp, (x) такую непрерывную 
функцию от x, что О<ы,(х) = 1 в В, в, (x)=), если x 
не принадлежит U, и в, (х) =1, если x принадлежит („.,; 
такие функции существуют, поскольку замыкания U,,,, ле- 
жат в U,. Тогда при всех x из (/, имеет место неравенство 


[ор ®) Фик (©) [< 272 [аифн (%) |. (ly 


В самом деле, это неравенство имеет место по построению, 
если х принадлежит И; если же х лежит вне Ч то 
левая часть неравенства обращается в нуль. Пусть 


со 


f(x) = Ханн, (х), ©). (2) 


1.7] СУММИРОВАНИЕ АСИМПТОТИЧЕСКИХ РЯДОВ 35 


В силу peice hi (1), этот ряд сходится для всех X и за- 
дает функцию f(x) в К (для x, не принадлежащих всем U,, 
этот ряд т Чтобы доказать, что fr Dann когда 
х—х., зафиксируем № и выберем x в пересечении Uy,, и 
R. Тогда [и (Х) ==1 при п==1,..., М и, в силу соотноше- 
ния (1), 


№ ©. © 
и— У а, | = У [арын | = | eva на | > ит 
N+1 Ат 


n=) 


= 2| Ones Физ: |= (py): 


Таким образом, )a,9, является асимптотическим разложе- 
нием до любого члена функции /, определенной равенством (2). 
Асимптотическая сумма ряда У, 4,9, является классом всех 
функций, асимптотически равных /. 

Окрестности (/, можно построить так, чтобы единствен- 
ной точкой, принадлежащей всем этим окрестностям, была 
точка X,; в этом случае ряды в формуле (2) обрываются при 
всех х-Ех,. Если функции $, непрерывны в К, то и функ- 
ция f непрерывна в R. Если х— вещественное переменное 
или точка П-мерного евклидова пространства, то функции в, (x) 
можно выбрать так, чтобы они были бесконечно дифференци- 
руемы; в этом случае, если функции Y, имеют непрерывные 
производные до А-го порядка включительно (Roo), тои 
функция f(x) имеет непрерывные производные тех же поряд- 
ков. Карлеман доказал, что для некоторых аналитических 
функций 9%, комплексного переменного x асимптотическая 
сумма содержит функции, являющиеся аналитическими функ- 
циями от х. 

В общем случае невозможно приписать асимптотическому 
ряду единственную асимптотическую сумму; но в некоторых 
случаях может оказаться, что при более точных предположе- 
ниях о коэффициентах асимптотических рядов и некоторых 
ограничениях на функции / (x) получаются единственные сум- 
мы; в этих случаях асимптотические ряды, несмотря на то, 
что они расходятся, в некотором смысле суммируемы к их 
асимптотическим суммам. Теоремы об асимптотических сте- 
пенных рядах, суммируемых к аналитическим функциям, ре- 
гулярным в некоторых секторах, были получены Ватсоном 
(1912а) и Неванлинна (1916). 
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ГЛАВА II 
ИНТЕГРАЛЫ 


Существует много методов для получения асимптотических 
разложений функций, заданных определенными интегралами. 
Копсон (1946) дал обзор этих методов; дальнейший материал 
содержится в лекциях ван дер Корпута и в статьях, указан- 
ных в конце этой главы. 


2.1. Интегрирование no частям 


Асимптотические разложения часто могут быть получены 
с помощью повторного интегрирования по частям, В качестве 
примера рассмотрим функцию f(x), заданную в области 
— п ао x < п интегралом 


м 


а я dt. (1) 
Повторно интегрируя по частям, о 
а 
а |p хх! (em arts. 
© > 
= Ох" (— 107 п рн.) 
n=0 о 


Можно доказать, что последний интеграл есть O(1), когда 
x— 0 в области Sy, А > —п/2. Мы получили, таким обра- 
зом, новый вывод эйлерова асимптотического разложения, 
рассмотренного во введении. 

Область применения этого метода довольно ограничена, 
и не всегда легко сформулировать точные теоремы, обладаю- 
щие достаточной общностью. В дальнейшем будут приведены 
некоторые результаты, которые нам кажутся основнымя. 
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Для любой функции f(t) обозначим через /„ производную 
т-го порядка, а через f_,, m-ii повторный интеграл. Таким 
образом, 


ат а 

=f Sun oi m==1,2, ..., (3) 
Чт __ a 

“Ff Sama а re eee (4) 


Заметим, что /_„ содержит M постоянных (по одной OT каж- 
дого интегрирования). Мы будем предполагать, что эти по- 
стоянные выбраны соответствующим образом. С помощью 
повторного интегрирования по частям получаем формулу 


8 N-1 
феи = У, s+ Ry © 
где г 
Sp = (— 1)" [&, (В) „-, (8) — в, (AR (MS (6) 
в 
R,=(— 0" | gy he, (Oat. 7) 


a 


Если (a, 8) является конечным интервалом, то формула (5) 
справедлива при условии, что функция g имеет непрерывную 
производную N-ro порядка, а функция # интегрируема. Если ин- 
тервал (a, В) бесконечен, то все интегралы должны быть схо- 
дящимися и должны существовать пределы g(t) Й_„_, (#} 
при {— 9, #— В. 

Если функция g имеет непрерывную производную М№М-- 1-го 
порядка, то с помощью еще одного интегрирования по частям 
получаем 

Ry=Sy+ Юм. (8) 


В некоторых случаях с помощью этого соотношения можно 
сравнить остаток Ry с первым отброшенным членом sy. 

Если функции g uh вещественны, а Функции В\й_уи 
В+. 4 м_, Знакопостоянны на интервале 9 < = Ви имеют 
на нем одинаковые знаки, то К у имеет тот же знак, что и 5х» 
причем |Ry|<|sy|. Для доказательства достаточно восполь- 
зоваться формулой (8) и заметить, что в нашем случае Ry 
и Кл., имеют противоположные знаки и, следовательно, Ry 
И Sy имеют одинаковые знаки, причем 


[В ==м|— [Юма || 
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Пусть функция g вещественна, функция |h_y_,| 803- 
растает и Е, бу. Знакопостоянны на интервале а < = 
и имеют на нем одинаковые знаки; либо пусть функция 8 
вещественна, |h_y_,| убывает, а функции Sy, @ ул. знако- 
постоянны на интервале A<t<B и имеют на нем раз- 
ные знаки; тогда | Ю| = 2|5%|. 

Мы докажем это утверждение в случае, когда [h_y_,| 
является возрастающей функцией от Ёи 2—0, Ly,, 0. 
При этих предположениях имеем 


| Км |= | м. (8) | (gy (8) — gy (2) < 
<|h_y_, (8) 5 м8) | —|A-y-1 @) 89) [< 
= [п м-а (В) & (В) —h_y_, (2) (9) |, 


и, следовательно, | Юл. |< |5л|. Из формулы (8) вытекает 
нужный нам результат. Если gy<0, Sy,,<0, то заменяем 
g на —g. В случае, когда функция |A_y_,| убывает, надо 
заменить Ё на — . 

В качестве приложения этих результатов рассмотрим 
функцию f(x), определяемую формулой (1). Если x >> 0, по- 
ложим 

g(x, деть big = 1" 


В этом случае 2„й_„20 при всех #20 и, следовательно, 
0=(— 1)” ®,=—(-— 1)” 5. Если х— комплексное число, 
то наши результаты неприменимы. Однако, если в форму- 
ле (1) заменить Ё на #/х и соответственно этому положить 
g=(1+7', Axo ехр (—#/х), то при комплексных зна- 
чениях х, таких, что Rex > 0, функции g,, и 8„.., знакопо- 
стоянны и имеют противоположные знаки, а, | является 
убывающей функцией от Ё при #20. Следовательно, в этом 
случае | R,,|<2]s,,|. (На самом деле и в этом случае легко 
доказать с помощью формулы (7), что | В,|=<|5.|.) 
Предположим теперь, что подынтегральная функция в фор- 
муле (5), а следовательно, также и функции Sy и Ry, зави- 
сят от переменной х. Если {5„} — асимптотическая последо- 
вательность и если, кроме того, с помощью одного ‘из ранее 
полученных результатов (или каким-либо иным методом) можно 
доказать, что Ау==О (54), то формула (5) приводит к асимп- 
тотическому разложению интеграла до N-ro члена. Напри- 
мер, в случае разложения (2) {х”} является асимптотической 
последовательностью при х —+0; выше было доказано, что 
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А, |= 2151 для любого N в области Rex > 0; следова- 
тельно, (2) является асимптотическим разложением функции 
J (x), определенным в области Sy, А > 0. (Фактически, во 
введении было доказано, что асимптотическое разложение 
справедливо в более широкой области Sy, А > —11/2; это 
также можно установить с помощью оценки последнего ин- 
теграла в формуле (2).) 

Последовательность {s,} часто оказывается асимптотиче- 
ской в случае, когда A(t) имеет вид 


h(t)=k (x?). 


Обозначим через А_„(и) т-й повторный интеграл по пере- 
менной и от R(u). Из формулы (5) вытекает 


8 N-1 
{ Ut) Get dt = У, (хи, (B) Rg (BX) — 


— Sn (2) Rens (x) Ry. (9) 


Далее, если функции №_„(и) ограничены, а функции 
8n (B) А_„_, (BX) — By (2) А-„_,(@х) не принимают значений, 
лежащих в некоторой окрестности нуля, и если, кроме того, 
Ry допускает приведенную свыше оценку, то (9) является 
асимптотическим разложением при х-—+с®, причем область 
изменения х определяется оценкой Ry. Кажущийся более 
общим случай h(t) = [xg (f)] может быть сведен к рассмотрен- 
ному выше с помощью разбиения интервала (х, В} на части, 
на каждой из которых функция ф(Й} монотонна, и введения 
на каждой из этих частей 9 (t) в качестве новой переменной. 
Чтобы применить один из описанных выше критериев для 
оценки Юл, надо знать поведение производных вида 


an eal 
Ч” |e (t)] * 


Если g WO имеют производные одного и того же или чере- 
дующихся знаков, эта информация может быть получена без 
непосредственных вычислений из результатов об абсолютно 
монотонных и вполне монотонных функциях (см., например, 
Уиддер, 1941, глава IV). Общие теоремы такого типа были 
получены ван дер Корпутом и Франклином (1951). Наиболее 
важным приложением этих методов является асимптотическое 
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разложение интегралов вида 
b b 
Ve yemat, {едва 


а 


Более общие результаты, касающиеся функций вида g(t) = 
= (1 — а) 5. (t), где 5, (В) — функции, имеющие непрерыв- 
ные производные, были получены ван дер Корпутом (1934). 


2.2. Интегралы Лапласа 


Интегралы вида 


f= fe *t 9 (t)dt = {9} (1) 


о 


называются интегралами Лапласа. Такие интегралы встре- 
чаются при решении дифференциальных уравнений с помощью 
определенных интегралов и во многих других вопросах. Не- 
собственный интеграл в формуле (1) понимается как 


lim { е (1 dt. 


Поэтому функцию o(f) мы всегда будем считать интегриру- 
емой на любом интервале 0 = {= Т, Т < oo. Функция ф на- 
зывается принадлежащей к классу L(x,), если интеграл 
в формуле (1) существует при х==х, в указанном выше 
смысле. Известно (Уиддер, 1941, глава П), что для любой 
функции $ из L(x,) существует ¥ {5}, причем {0} является 
аналитической функцией oT x в полуплоскости Re (x — x,) > 0. 
В частности, если функция y(t) интегрируема на любом ин- 
тервале 0 == 7, Т« сю, и если ¢(t)==O(e?* ') для неко- 
торого постоянного a, когда #— с, то ®{9} существует 
(причем несобственный интеграл абсолютно сходится) и яв- 
ляется аналитической функцией от x в полуплоскости 
Re (x — a) >0. 

В некоторых случаях асимптотическое поведение f(x) при 
х—+ со может быть исследовано с помощью интегрирования 
по частям. Если g(t) имеет при 09=Ё=а непрерывную 
производную N-20 порядка и принадлежит, при некотором 
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x,, к L(x,), mo 
F(x) У" (Oxo (2) 
до М№-го члена равномерно по atg xX, когда X—»+ 00 6 Ss, A>0. 


В самом деле, пусть Ве (х —х,) > 0. Запишем интеграл (1) 


в виде 
а 


бе = Jet of att Se y(t) at. (3) 


cy 


Второй интеграл существует и приводится с помощью инте- 
грирования по частям к виду 

wo 

3 

{ -®-=9 1 (1) dt, 


a 


| ee 


a hy 


где функция 
t 


b(t) =} em то (и) аи 
а 


ограничена при ta. Если |%| < А, то 
Aem = 20) а 


lx—x ll? — | 


[1 [= == O(e~"), 


когда x ==9-+ © —* co B Ss, A> 0. К первому интегралу при- 
меним формулу 2.1(5), положив g= y(t), й_„==(— дея: 
Получим 


5,= [ye (0) —y (a)7*4] x7" ey (0) х-"-* +O (е-”) 
и 


а 
Вых [gi (бе dt. 
о 
Функция y(t) непрерывна и, следовательно, ограничена при 
O<t<a, Если |9% (В | < В, то 
|Ry| << B|x|-%e7* < В|х|7^-* совес А =0О(х7 474) 


равномерно по ато х, когда х—* со в 54, А`>0. Чтобы за- 
вершить доказательство, заметим, что O(e7**)==0(x~A) 
равномерно по arg x, когда х— co в Sa, ADO. 
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Существенное обобщение последней теоремы основывается 
на следующей лемме. 

Лемма. Пусть o(t) и Y(t) при некотором X, принадле- 
жат L (x,), b(t) > 0, = SHU {5}. Если eg (9) — со, 
когда р—>со при каждом а`>0, и если + =0($ (1), 
когда t—+0, mo f(x) =0(8(5)) равномерно по atg x, когда 
x==ptis—oos Ss, AD>0. 

Доказательство. Задалим ¢>>0. Найдется a> 0 
такое, что |2 |= = при O<Ct<a. Запишем & {$} в виде (3), 
где а имеет выбранное значение. Каки в доказательстве пре- 
дыдущей теоремы, второй интеграл есть Oe”) при p—oo и 

а. а 
[емо а | ce | ey (at egy) 
0 о 
Таким образом, 
If (*)| —. fe"? 
zy ee Oras? 
Ho при достаточно больших 6 правая часть этого неравенства 
не превосходит 2s. Поэтому /(х)==0(9(7)}. Равномерность 
по arg x вытекает из того, что в области S, имеем: 
|x [12 == созес А. 

Из этой леммы вытекает следующая теорема. 

Пусть при п=1, ..., N функции ,,(t) принадлежат 
при некотором x, к L(x, 4, (t) > 0 при i> Ou g,=L{9,}. 
Если {'»,\ является iy ia se последовательностью 
при t—+0 и е"'2, (2) —+ со при всех а`>0 и всех п, когда 
р —> oo, mo {2„(2)} является асимптотической последова- 
тельностью при в —+- со. Если при тех же предположе- 
ниях %(t) repaint dwesicah L(x,) и 


vo ~ da,d,() до М-го члена при 10, 


mo 
1(2) > a,g, (9) до N-ro члена при p—-} oo. 


Если, кроме того, при всех п==1,..., М функции g,(9)/g, 
ограничены в Ss при достаточно больших |x|, mo {g,(x)} 
также является асимптотической последовательностью и 
f(x) ~ Уав,(х) д0 N-20 члена равномерно по atg x, когда 
x— oo 6 Ss, А 0. 

Доказательство. Из леммы вытекает, что &„., (2) = 
=0(g,,(2)), и, следовательно, {Sn (6)} является асимптотической 
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последовательностью при 2 —> oo. Чтобы доказать, что разло- 
жение функции Y является асимитотическим при р—+с, 
заменим в лемме Ф на by и ф на 


т > 

— 5} 

9—2, 
a= 


Когда x—>co в Ss, TO, в силу леммы и дополнительных 
предиоложений о g, (x), имеем 


Ен: ов, )) = ов, (и) = 0 (в, (9). 


Следовательно, {g,,(x)} является асимптотической последова- 
тельностью при х—+ 00 B Sy. Дальнейшее доказательсгво 
протекает совершенно так же, как и выше. 

Наиболее важным частным случаем является 


b= tat, <<... 
Все условия теоремы выполнены; в частности, 
8 (х) =Г(,) хи 
и потому для всех x из 5,, A>>O, выполняется неравенство 
< (1) овес 


Мы получаем, таким образом, следующую теорему об асимп- 
тотических рядах степеней. 

Пусть (< <A, <<... « Если g(t) при некотором x, 
принадлежит L(x,) и 

g~ Уа,Ё"-: до N-ro члена при #—0, 
то 
/^ УГО) а,х-» до М-го члена равномерно по ах 
при x—oo в Sys, ADO. 


Отметим еще следующие примеры асимптотических после- 
довательностей, к которым применима общая теорема; 


ws ents __ (a —1)! 

b,=(1—e "3, Е. (и И’ a 
dni. Gui не =) 

b, == (e — 1)" -*, 8u= и. “по, г 


lo, Е" Qn — 2)! 
%=(2% 5) » 8 Goa pnw Fy. panty 4) 
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В случае асимптотических степенных рядов, равно как и 
в случае (4), N может быть как конечным, так и бесконеч- 
ным; в случаях (5) и (6) N должно быть конечным. ь 

Результат, полученный с помощью интегрирования по 
частям, является частным случаем теоремы об асимптотиче- 
ских рядах степеней. Если функция ¢(f) имеет при 0 == 
непрерывные производные до N-ro порядка включительно, то 
из формулы Тейлора с остаточным членом в форме Пеано 
вытекает, что при #—+0 

о ~ So (0) "jn! 
до N-ro члена. Поэтому (2) вытекает из теоремы об асимп- 
тотических степенных рядах. 

Во многих случаях можно расширить область, в которой 
справедливо асимптотическое разложение, до области вида Sy, 
А< 0. Если $(#) — аналитическая функция от № регулярная 
в Si, и если y(t) == O(e* Hy для некоторого @, когда #—+ со 
B Ss, то можно повернуть путь интегрирования в формуле (1) 
так, "чтобы он стал лучом, лежащим в 5, Таким образом, 
функция F(x) аналитически продолжается в ‘некоторую область, 
которая содержит ceKkTop —п--6<Са!в (х —а) “п—0 
(см., например, Дёч, 1950, стр. 362 и далее). Если {%,} 
имеет соответствующие свойства вдоль каждого луча и 
$ Dadar когда # —0 в $,, то f~ Уа,<.„, когда x — co 
в 54, где 4 >>4—n/2. Сформулируем точную теорему для 
асимптотических степенных рядов. 

Пусть $ (8) — регулярная функция в S,, y= O0(e") pas- 
номерно no argt при некотором а и #00 в S,, 


ф-— da,t-* до М-го члена 
равномерно по argt, когда #—08 S,, 
<, <}, = вы <hy 
Тогда f(x) существует в секторе 


—n-+6< arg (х — а) п— 0 


~ Уа,Г(,)х-”" до М-го члена 
равномерно по агвх, когда x—> со в Ss, AD>A—n/2. 


где 


Частный случай этого результата, когда числа 2, обра- 
зуют арифметическую прогрессию и ф может быть представ- 
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лено в виде сходящегося ряда 


Ура, -1 

для достаточно малых |Ё| в S,, известен как senna Ват- 
сона; эта лемма достаточна для большинства приложений. 

На протяжении этого пункта было изучено поведение 
интеграла (1) при болыних значениях x. Подобные методы 
могут быть использованы для исследования функции f(x) 
при х—>х,. Относительно основной леммы см. Эрдейи (1947), 
а относительно некоторых наиболее важных результатов 
см. Дёч (1950, глава 13) и Уиддер (1941, глава 5). 


2.3. Критические точки 


Мы видели в последних двух пунктах, что при некото- 
рых условиях асимптотическое поведение интегралов опреде- 
ляется поведением подынтегральных функций в определенных 
точках; например, в рассмотренных выше случаях такими 
точками были концы интервала. Подобные точки ван дер Кор- 
пут (1948) назвал критическими. Мы изложим сейчас некото- 
рые методы, позволяющие находить асимптотические разло- 
жения, соответствующие основным типам критических точек. 

Рассмотрим сначала интегралы вида 


8 
еде at, (1) 


a 
где x — большой положительный параметр и A(t) — вещест- 
венная функция. Если h(t) имеет максимум в точке f==t и 
h(t)<Ch(t) при t=4%, то при больших значениях X модуль 
подынтегральной функции имеет резко выраженный максимум 
в точке, весьма близкой к точке т. Поэтому поведение ин- 
теграла (1) при больших значениях х зависит в основном 
от значений функций g(t) и A(t) в окрестности точки <. 
Интеграл может быть приближенно вычислен с помощью 
разложения этих функций в окрестности точки Tt. Это является 
центральной идеей метода Лапласа (см. п. 2.4). Мы столк- 
нулись с подобным случаем в п. 2.2, где #(=—Ь 
0=1< с, и потому A(t) имело максимум при #=0. 
В соответствии с этим мы нашли асимптотическое разложение 
интегралов Лапласа путем разложения функции g(t) при ма- 
лых значениях &. 
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Если x и A(t) комплексны, причем g(t) и A(t) являются 
аналитическими функциями от Ё, то часто оказывается воз- 
можным изменить путь интегрирования так, чтобы он прошел 
через одну или несколько точек, в которых h’ (1) =0. Если 
<—-одна из таких точек, то она является критической; 
в этом случае можно провести часть пути интегрирования, 
лежащую в окрестности точки т так, чтобы x [A (В —# (®)] 
было вещественным на этой части. После этого интеграл 
может быть вычислен с помощью метода Лапласа. Опи- 
санный метод называется римановым методом перевала 
(см. п. 2.5). 

Перейдем теперь к интегралам вида 


3 
\ еде at. (2) 


Предположим снова, что х-— большой положительный 
параметр, а А(!) — вещественная функция. В этом случае 
из-за быстрых колебаний функции exp [1хА (1)] происходит 
интерференция значений подынтегральной функции. Это не 
имеет места лишь на концах интервала и в стационарных 
точках функции A(t). Поэтому в данном случае критическими 
точками являются концы интервала’ (%, В) и стационарные 
точки функции #(1)*). Если k(f) не имеет стационарных 
точек на интервале << В, то хорошее приближение, как 
правило, получаегся с помощью интегрирования по частям 
(п. 2.1). Для оценки вклада, получаемого от окрестности 
стационарной точки T, применяется метод стационарной фазы 
Стокса (п. 2.9), основанный на разложении функций g uh 
в окрестности этой точки. 

Метод стационарной фазы был распространен ван дер 
Корпутом (1936) на интегралы вида (1), где xh(t) может 
быть любой комплексной функцией (в методе стационарной 
фазы эта функция должна быть чисто мнимой). Согласно 
ван дер Корпуту, в этом случае критическими являются 
точки т, в которых функция x/2h' (1) [1" (1]]-'з веществен- 
на, причем мнимая часть этой функции меняет знак при 


*) Легко видеть, что около концов интервала не происходит 
полной интерференции значений подынтегральной функции, а в ста- 
ционарных точках функции A(t) функция ехр [1х (#)]| медленно ме- 
няется. (Прим. перев.) 
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переходе # через +. (В случае интеграла (2) h=-ik, где 
k— вещественно, а потому единственными точками, в кото- 
рых функция 

дв’ (ву == хи" А") 


вещественна, являются стационарные точки функции А.) 


2.4. Метод Лапласа 


Пусть в интеграле 
3 


f(x) == ( вета (1) 


a 


(К) является вещественной функцией вещественного пере- 
менного #, функция &() может быть как вещественной, так 
и комплексной, а х — большое положительное переменное. 
По Лапласу, основной вклад в интеграл дают окрестности 
тех точек интервала х=2=<В, в которых A(t) принимает 
наибольшее значение. Если A(f) имеет конечное число макси- 
мумов, то интервал (я, В) можно разбить на конечное число 
интервалов, на каждом из которых функция (tf) принимает 
максимальное значение в одной из концевых точек и не при- 
нимает его в других точках. Поэтому мы можем считать, 
что в интеграле (1) функция # (1) имеет максимум при =, 
причем A (2) < h(a) при я < Е= 8. 

Предположим, что функция g непрерывна, а функция й 
обладает непрерывной производной второго порядка, h’ (а} = 0, 
h" (a) < 0. Введем новое переменное и, определяемое равен- 
ством й(5)— (В ==и?. При достаточно малом 1 функ- 
ция A’ (2) отрицательна в интервале х < 1< я -[- 1. При х-— оо 
имеем 

ат 
| 

fr } ве" dt= 
a 


=—| ou 2 texp x [й (2) — и*]} du, 
0 


где (= [й («) —# (= -|-1)|'*>0. Так как существенным яв- 
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ляется лишь значение подынтегральной функции в окрестности 
точки “=O, мы можем приближенно заменить g(t) на g(a), 
а выражение и/й’ (#) — его пределом при f—+a, то есть зна- 
чением — [-— 2h" (2)|-"*. Мы получим при этом 


—2 


109 — [ml "в oa) ( {ехр [— хи? + xh (a)]} due. 


a6 


В силу тех же соображений, можно заменить интервал интег- 
рирования интервалом [0, со). Это приводит к формуле 
Лапласа 


fore ayer |e 


(2) 
Буркхарлт (1914) и Перрон (1917) показали, что этот pe- 
зультат может быть доказан с помощью разложений функций 
хи йв окрестности точки ©. Копсон (1946) воспроизвел 
простое доказательство Пойа и Cere, a Уиддер (1941, 
глава VII) дал более сложное доказательство при более общих 
предположениях. Дальнейшие обобщения формулы Лапласа 
были получены Су (1949а, b; 1951а, b), Леви (1946) и 
Руни (1953). Метод Лапласа был использован для интегра- 
лов, зависящих от двух больших переменных, Фальксом (1951) 
и Томсоном (1954), а к двойным и кратным интегралам при- 
менен Су (1948a, b; 1951с) и Руни (1953). 

Приводимое ниже обобщение результата Лапласа может 
быть выведено из установленных выше свойств интегралов 
Лапласа. Пусть g и h— функции, заданные на интервале 
(a, В), такие, что интеграл (1) существует при доста- 
точно больших положительных значениях х. Пусть, далее, 
функция h вещественна, непрерывна при Ё==а, имеет не- 
прерывную производную при в <= а- 1, 4>0, причем 
WO при a<t<a+y, h(\)<h(z)—e, #20, при 
&--1=2= В; наконец, пусть h'(t)~-—a(t—a)-' u 
g(t)~o(t-— <)" при t—+a, где 4>>0, 9>0. Тогда 


8 
Пе == exit ai~+T (+) (=) et, (3) 
$ ¢ 
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Для доказательства заметим сначала, что 


| | вое |= 
а} 


8 
<exp {x [h(a)—e]} | eat 
a+y 
= 0х №2], х-ьс. (4) 
На интервале (a, a = 1) введем новую переменную 
и=# (2) —#( и положим И=й (а) —1№(&-1)>0, 
& (и) = — g (d/h (0). Мы получим 
aty U 
Y g(t et at =o | в (иде ди. (5) 
é 


@ 


Ho 
t 


и (а) В) ож S¢—a), ta, 


и потому 
imax (")". и—0. 
Итак, 
а, toa, 
и, следовательно, 
kw ~2 (2), u—0. (6) 


В силу полученных выше результатов об асимптотическом 
поведении интегралов Лапласа, из соотношений (5) и (6) имеем 
«1 
\ g (tye! dt~ 


By 
В Гу \№ h = 
-у (= ) Г (> ) хе" ®, ¢ + 00. (7) 


Из соотношений (4) и (7) и вытекает асимптотическое равен- 
ство (3). Заметим, кроме того, что соотношения (4) и (7), 
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а следовательно, и соотношение (3), остаются справедливыми, 
если х является комплексным переменным и x—>0o в Sy, 


A>0. 


Покажем теперь, как получается асимптотическое разло- 
жение функции f(x). В нижеследующих формулах положено 
п=0,1, ..., М— 1. Если 


—h' )~ Уд)", 
g(t) ~ Yb, (tay 


до N-ro члена при #—+&, то справедливо разложение 


(8) 


g(t) Leven 
—FH ^ Dealt—%) 2 до N-ro члена при t—+a. (9) 
Коэффициенты c, в этом равенстве могут быть вычислены 
с помощью формального деления. Таким образом, 


1 
и=— (м @a~y tay" 
до N-ro члена при t—+a., (10) 


Из полученного выражения вытекает, что “YY может быть 
разложено в асимптотический ряд по степеням #— а. Под- 
ставляя этот асимптотический ряд в соотношение (9), полу- 
чим асимптотическое разложение вида 


& (и) — Уи" -» до М-го члена при и—+ oo. (11) 


Используя соотношение (11) вместо соотношения (6), по- 
лучим асимптотическое разложение до N-ro члена для инте- 
грала (5); кроме того, оценка (4) может быть заменена более 
точной оценкой о (x7 O+M/%e** ), Таким образом, 


До ~ e@@ У уг (+27) 2-o-+m (12) 


У 


до N-ro члена, когда х—+ ов Sy, A> 0. Коэффициенты y, 
могут быть вычислены с помощью формальной подстановки 
согласно описанной выше схеме. 

Существует другой метод вычисления коэффициентов y,, 
позволяющий избежать обращения асимптотического ряда (10) 
для получения разложения и" по степеням #—%. Именно, 
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из соотношения (10) имеем; 


h(t) = h(a) — ад" №, (0, 


где 
N-1 
an эт v4N- 
h(t) = — SY ев)" Ноев). 
ast 
Полагая 
a+ aby 


| g(t) er at ой | I(t) exp [9—5] dt, 


а 


разлагая формально функцию 
1 (t)= g(t) exp [xh, (4)] 


по степеням ¢#— и интегрируя почленно, мы приходим 
к разложению (12). Можно доказать соотношение (12), ис- 
пользуя эту схему. 


2.5. Метод перевала 


Рассмотрим интеграл 
8 
Ду у gihe™ at, (1) 


где x — большое комплексное переменное, © и A — аналити- 
ческие функции комплексного переменного # а интеграл 
берется вдоль некоторого пути в комплексной Ё-плоскости. 
Этот интеграл может быть асимптотически вычислен с по- 
мощью метода перевала, впервые введенного Риманом и раз- 
витого Дебаем. Копсон (1946) дал подробное описание этого 
метода, приведя многочисленные ссылки и примеры. 

Точки #-плоскости, в которых A’ (1 ==0, называются mo4- 
ками перевала или седловинами. Поверхность, изображаю- 
щая jexp[xh(t)] как функцию от Ё и 1щЬ называется 
рельефом функции е”"; на этой поверхности точки, в кото- 
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рых A’ (1 =0, являются точками перевала, и наиболее удоб- 
ный путь из одной «долины» в другую ведет через один 
или несколько перевалов. Пусть t является точкой перевала; 
если fh’ (т) =” (=... ==” (1) =0 и AME (5)520, то 
< называется точкой перевала порядка т. Кривые на пло- 
скости ¢t, вдоль которых Ве хй (ft) имеет постоянное значение, 
называются линиями уровня: вдоль таких кривых функция 
е^" имеет постоянный модуль (эти линии являются горизон- 
талями рельефа), фаза же функции e** меняется вдоль линий 
уровня быстрее всего. Линии, вдоль которых Imxh(t) по- 
стоянио, называются линиями стока: вдоль таких линий е* 
имеет постоянную фазу, а jer” | изменяется быстрее всего 
(эти линии являются градиентными линиями рельефа). В точке 
перевала х порядка m пересекаются под равными углами 
т-{-1 линия уровня; эти углы делятся пополам ш-|-1 ли- 
нией стока; вдоль каждой из линий стока т является ста- 
ционарной точкой функции |'е** © |. 

Метод перевала состоит в том, что путь интегрирования 
деформируется так, чтобы он, насколько это возможно, со- 
стоял из дуг линий стока. Если % и В лежат на линии стока, 
проходящей через перевал, например, если а и В являются 
особыми точками функции A(t), то путь интегрирования 
можно деформировать так, чтобы он полностью совпал 
с линией стока, проходящей через перевалы; в противном 
случае могут встретиться две линии стока, не проходящие 
через перевалы. Последний случай может быть описан с по- 
мощью рельефа следующим образом; сначала мы спускаемся 
вдоль линии стока в особую точку, а потом вдоль другой 
линии стока поднимаемся на перевал. Функция xh (1) моно- 
тонна вдоль всей линии стока (исключая перевалы), а потому 
для асимптотического вычисления интеграла может быть 
применен метод Лапласа. Асимптотические разложения функ- 
ций g и hf, необходимые для применения теоремы предыду- 
щего пункта, суть не что иное, как тейлоровские разложе- 
ния этих функций в окрестности точки линии стока, в которой 
хй(Ё) достигает наибольшего значения (эта точка часто 
совпадает с точкой перевала). Обращение ряда 2.4(10) может 
быть выполнено с помощью формулы Лагранжа (см., напри- 
мер, Маркушевич, 1950, стр. 342). 

Мейер (1933a, b) показал, что оценка погрешности мо- 
жет быть получена с помощью разложения по формуле 
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Лагранжа с остаточным членом. Он показал также, что 
в некоторых случаях можно получить рекуррентные соотно- 
шения для коэффициентов. 

Мы рассмотрим сейчас ряд примеров на применение ме- 
тода перевала; эти примеры заимствованы из работы Коп- 
сона (1946). 


2.6. Интеграл Эйри 


Мы изучим сейчас асимптотическое поведение интеграла 
wo 
Ai(e)==- | cos 4 st-Les\ ds (1) 
® 3 р. 
0 


при больших значениях 2. Подстановка 


2112$, (2) 
приводит этот интеграл к виду 
Я xe $ в. al 
Ai(x**) = 2 exp [ix (g#+#)| dt. (3) 
Со 4 


К этому интегралу может быть применен метод перевала. 
Будем рассматривать в интеграле (3) # как комплексное пе- 
ременное интегрирования. При х`>0 путем интегрирования 
является вещественная ось #. Этот путь можно деформировать 
в кривую, которая начинается в бесконечно удаленной точке, 
проходит по сектору 21/3 < а=е{Ё< п, переходит в сектор 
O<arg¢<(m/3 и заканчивается в бесконечно удаленной 
точке. В нашем случае 


в0=: (зе), h(t) =i +1). 


Поэтому перевалами являются точки #==-- 7 (нули функции 
в’ (1)). Линии стока определяются равенством Im A (f) = const. 
Полагая #=&-|-Й, получим 


ши =З — Of +6, ши (д =0. 
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Поэтому уравнение линии стока имеет вид 
Е — 377 +3) =0. (4) 


Это уравнение вырожденной кривой третьего порядка, со- 
стоящей из мнимой оси и двух ветвей гиперболы. На рис. 1 
стрелками обозначены направления, вдоль которых Re A(t) убы- 
вает. Асимитотами гиперболы 
являются линии §+4V3=0, 
и, ‘очевидно, путь  интегри- 
рования в выражении (3) мож- 
но деформировать так, чтобы 
он превратился в верхнюю 
ветвь гиперболы и wea, 
таким образом, из точки 
со -ехр (5116) в oo-exp (im/6). 
Вдоль этого пути интеграл (3) 
сходится, если Rex > 0. 
Положим 
ЭДА Рис. 1. 
оо+ехр (ir/6) ор-ехр (5й*(6) 
= J — J a—1—1, 6) 


é é 


и применим к интегралам J,, /, метод Лапласа. В обоих 
интегралах функция A(?)—A(i) вещественна и достигает 
максимума в точке #== следовательно, #(А) —й(1) яв- 
ляется убывающей функцией. Введем новую переменную и, 
положив 


ив) — п = — 2 (зе) = 
ЕВ С 


Из (6) вытекает 
Чи i [1 eit д |, (7) 


где и — положительный квадратный корень, а [ в jp? — 
ветвь функции, которая равна единице в точке 2—4 
верхний знак в (7) относится к J,, а нижний знак — к /,. 
Из теоремы Лагранжа вытекает, что в достаточно ма- 
лой окрестности точки перевала ¢—i обладает разложением 
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вида ¢—i== Sb, (5 ")", где пб, — коэффициент при 


{{—1)"-' в разложении [1 —i(¢—/3]-" по степеням 
1—1. Таким путем получаются разложения 


(8) 


где верхние и нижние знаки относятся, соответственно, к /, 
и /,. Ho 


нм 
= [ee du 
о 


и, в силу п. 2.4, асимптотические разложения интегралов 
/, н Г, получаются путем подстановки разложений (8) в 
а|4и и последующего почленного интегрирования. Таким 
образом, 

© 


У (= 1)" oP (3n/2 — 1) 


ao 
exis], = ae oa ge Ц"! du ~ 
ь 5 


fora} 


= 10" P-L (8n/2 — 1) 


is 2 — 1131 


п 


Подставляя это разложение в (5) и выражая полученный 
результат через 2, получаем после некоторых упрошений 
разложение 
@ 
1 OFX Pm + 1/2) ae 
; РИ: 29 —— — g3/2 I (— 923/)-т 
Аа ~ ex ( 2 p> Om (— 929/)"™, (9) 
то 
Это асимптотическое разложение справедливо равномерно 
по ага 2, когда Z—+oo в области [ага 2| < п/3 —А, А>0. 


2.7. Дальнейшие примеры 


Мы рассмотрим теперь два примера, в которых пределы 
интегрирования не являются особыми точками и поэтому 
асимптотические разложения не могут быть получены с по- 
мощью разложений в окрестности перевала. Кроме того, 80 
втором примере перевалы имеют второй порядок, 
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Первый пример таков. Пусть х>0 и 
ao 
(x) = | exp [ex (s#+4)| dt. (1) 
о 


Функция A(t) в этом примере та же самая, что и в п. 2.6, 
а потому линии стока совпадают с изображенными на ри- 
сунке в п. 2.6. Легко видеть, что соответствующий путь 
интегрирования состоит из части мнимой оси между точками 
О иги половины верхней ветви гиперболы. Таким образом, 


é oo-exp (ins) 


= [ etme, (2) 


i 


Асимптотическое разложение второго интеграла было уже 
получено выше. В первом интеграле функция A(t) вещест- 
венна и убывает, когда Ё изменяется от 0 до i, поэтому 
снова можно применить метод Лапласа. В соответствии 
с этим методом положим 


и—1(0) —1 (8 —# (1 sf). 


3 
Из теоремы Лагранжа вытекает, что —й= Ур", где 
nb, является коэффициентом при (—Й)”`* в разложении 
(1 ив |3)-" по стеиеням — Й. Очевидно, что b,==0, если 
п четно и 


{Зт т +1 
2 т! (2m -- 1)13" * 


Подставляя это выражение в первый интеграл равенства (2) 
и интегрируя почленно, получаем 
т 
ее dtm 
e 
ao fos) 
Зт)! и?" P | 
xa yy бт. du ~i Sy т 


m! (2m)! 3” т! 3” * 
mao m=0 


Ц 
®— 


Из результатов п. 2.6 вытекает, что второй интеграл фор- 
мулы (2) экспоненциально мал по сравнению с первым. Мы 
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получили, таким образом, следующий результат: 


oO ao 
{exp [ix (+ +4)! dt~i > eee (3) 
о 


т —=0 
при х—+ сю в Ss, ASO. у 


В заключение рассмотрим интеграл 


F(x) == § exp (ixt*) de, 4) 


0 


tae х>0. Здесь A(t) it? и #==0 является перевалом 
второго порядка. Линии стока, проходящие через перевал, 


Рис. 2. 


задаются равенством Im (#*) ==0, то есть являются линиями 
t=n/6, Е т/2, 51/6. На рис. 2 стрелками указаны 
направления убывания функции |exp(éxf*)| при x >0. Ни 
одна из этих линий стока не проходит через точку #==1. 
Уравнение линии стока, проходящей через точку # ; 
имеет вид Im (éf*)==1. Полагая #==&--Й, можно записать 
это уравнение в виде & — 31° ==1. Это — кривая третьего 
порядка, ветвь которой, проходящая через точку #=1, 
изображена на рисунке. Для того чтобы перейти из точки 0 
в точку 1 вдоль линий стока, мы интегрируем сначала от 0 
до со вдоль линии агоЁ==п|6, а потом от со до 1 вдоль 
верхней половины ветви кривой третьего порядка. В соот- 
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ветствии с этим мы получим 


ор-ехр (xi) 6) w-exp (nf6) 


f= J — J срывал, 6) 
о 1 
BJ, положим z==— if? или == e's, где uv? >> 0, и 
получим 
fea) 
J а ge з4и==Г (4/3) ©1716 х- 1, (6) 
о 
В J, положим и==—2(Ё — 1) или #==(1-- Ш) и получим 


fo} 
зе жи (1 —|- м) 2 ди. 
о 
Разлагая (1--Й-?? no формуле бинома, имеем 


I~ rms el* Утеча 3) (ix)? (7) 


Подставляя (6) и (7) в (5), получаем, наконец, 


1 


§ exp (де) dé ~ I (4]3) в x -¥8 — 


— ro =". Ги 23) ()-"-* (8) 


п=о 


при x —+co в 54, А 0. 

Последнее равенство описывает асимптотическое поведе- 
ние функции f(x), когда х—» co в правой полуплоскости. 
Если x —> + joo, то подынтегральная функция в равенстве (4) 
вещественна и можно применить метод Лапласа. Наконец, 
если х—» со в левой полуплоскости, то можно применить 
соотношение 


I (x)=f(— x), 


вытекающее из (4), где черточка означает переход к комп» 
лексно сопряженному выражению. 
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2.8. Интегралы Фурье 


Интегралы вида 


еда (1) 


2 
ta 


называются интегралали Фурье. Мы будем считать здесь, 
что (4, 8) является вещественным интервалом, как правило, 
конечным; функция $Ф(!) предполагается интегрируемой, так 
что интеграл (1) существует при всех вещественных х. Мы 
изучим асимптотическое поведение интеграла (1) при x —-[- oo; 
чтобы получить асимптотическое поведение при x —+ — co, 
достаточно заменить # на —. В отличие от интегралов 
Лапласа (п. 2.2), в рассматриваемом случае единственным 
эффективным методом для получения разложений является 
повторное интегрирование по частям; лишь в случае, когда 
y(t) — аналитическая функция, может быть применен метод 
перевала. 

Докажем сначала следующее утверждение. Если функция 
y(t) имеет на отрезке а << {=5 В непрерывные производные 
до N-20 порядка включительно, то 


в 
{ ely () di = Ву) — Ао (хм), x00, (2) 


20e 
N= | 
Ay(x) = > GO (gh Bel, 
n=0 
Naa | (3) 
By (x) mes ee (В) х-"-1е 
= 
ay™==d"yldt". Этот результат остается справедливым, 
если &=— со (или В==оо) при условии, что 4 (1—0 
при t—+— co (или t—+ 0c) для всех п==0, 1,..., М— 1, 
причем oN’ (1) является интегрируемой функцией на (<, В). 
Чтобы доказать это утверждение, применим формул 
2.1(5), положив ==, g,=9", ho=e™, и „== (1х)-" ей. 
Для остаточного члена получаем 
Ry==(— ix | ем (у. 
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Так как при х—* со интеграл по лемме Римана стремится 
к нулю, то Ry=o0 (x7), 

Заметим, что если функция ф и ее первые М— 1 произ- 
водных обращаются в нуль в точке а (например, если Ф 
тождественно обращается в нуль в некоторой окрестности 
точки a), TO Ах(х)=0; точно так же Ву (х) = 0, если функ- 
ция фи ее первые М— 1 производных обращаются в нуль 
в точке В (например, если ф обращается в нуль в некоторой 
окрестности точки 5). 

Перейдем теперь к рассмотрению интегралов Фурье, 
в которых подынтегральная функция имеет особенности 
простого типа на одном из концов интервала. 

Пусть функция g(t) имеет на отрезке и < &= В не- 
прерывные производные до М№-го порядка включительно; 
y™ (8) =0 при п=0,1, ..., N—1 и 01. Тогда 

4 

ето — Ay (x) 0-м, x -+ 00, (4) 

2 
20e 

М1, ы 
дну У, EE ое ах. (5) 


n=0 


Пусть функция ot) имеет на отрезке o<t< В не- 
прерывные производные до М-го порядка включительно; 
gy (2) =0 при п==0,1, ..., М—1и O< p<]. Тогда 


ев 2)" офа = By(x) +O"), x-> 00, (6) 


а 


где 


N-1 
By (x)= У ПЕН в) ert (п — вр (В) х-"- Рег. (7) 


в=о 


При \==1 равенство (5) превращается в первое из pa- 
венств (3), а при и==1 равенство (7) превращается во вто- 
рое из равенств (3); однако О-члены в (4) и (6) дают мень- 
шую информацию, чем о-члены в равенстве (2). Вместо 
О (х-^) можно писать o(x7N7**") в (4) и o(x4-**) в (6); 
эта форма записи остастся справедливой при k= 1 или в =1 
соответственно. 
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Докажем равенство (4). Применим формулу 2.1(5), поло- 
жив «(д = (0, £, =H), Al) =A, (ев) и 


ВА (шим du, (8) 
1 
n==0, 1, ..., N—1. 


В равенстве (8) мы считаем Ё`> я и выбираем путь интегри- 
рования, целиком расположенный в квадранте 0 <arg (и — 2) = 
< п|2. Интеграл абсолютно сходится, и 
а 
ai tan-1 (= „(0 по %... Nt 
Если выбрать в качестве пути интегрирования луч и ==#-[- is, 
o>0, то получим [и—а| =, |u—af?<(t—a)? 
при 0<\-=1; следовательно, 
+ ю5 
\ ju—t|"[e** dul. 
` 
1 


1 
(t(~ay? 
п! 


hy. |< 


Подставляя u==t-+-is, получаем 
jh, |< (t—a)P tx"; ра, х>0. (9) 
Далее, из равенства (8) вытекает 


(— пр afin 


h_y-, (@) = | (u—a)"**-? e!** du, 


п! 
и, полагая u==a-l-is, >20, имеем 

ое 1 ED очен", (10) 
n=0, 1, see 


Мы можем теперь применить формулу 2.1(5). Члены в $» 
зависящие от В, обращаются в нуль, так как 4” (В) =0. 
Из (10) вытекает, что \5„=-—Ал, где Ay задается фор- 
мулой (5). Кроме того, из 2.1(7) следует 


Ry=(- 0% (ВВ y(t) dt, 


¢ 
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и потому, в силу (9), 
8 
[Ryl<a-¥ | [и (t— ay" dt = Oe), 
а 
Тем самым равенство (4) доказано; равенство (6) доказывается 
аналогично. 

Рассмотрим, наконец, интегралы Фурье, в которых подын- 
тегральные функции имеют особенности на обоих концах 
интервала. 

Если функция g(t) имеет на интервале я < #=В непре- 
рывные производные до N-ro порядка включительно и 
Oil, O<p<l, то 


8 
ив) в би 


== By (x) — Ay(x)-+ О(х-^), х-ьо, (11) 
где 
Ay (x)= 


N-1 & 
= У i es ) e В Ee [8 — «2-1 + (2)], (12) 
по 


N-1 
By (x)= У Tete) ent eo— ag ть [8 — 8) ф (8). 
n=0 


Если h=p=1, то О(х-^) в формуле (11) можно заменить 
на о(х-^). 

Эта теорема содержит три полученных ранее результата 
как частные случаи. Чтобы доказать равенство (11), мы 
используем часто применяемый для аналогичных целей прием. 
Введем функцию y(t), которую ван дер Корпут называет 
нейтрализатором. Эта функция y(t) бесконечно дифферен- 
цируема на интервале & = { == В, причем у (0) = 1, У” (a) =0, 
n=1, 2, ..., У? (8) =0, п==0, 1, 2, ... Примером такой 
функции является 
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Пусть ¥(f) — нейтрализатор. Тогда 


8 
миа (B— об 


8 
= J el ¢@— a) fy — да 


++ § ig — Af" о 91] ay pip at. (18) 


Первый интеграл в правой части этого равенства получается 
из интеграла (4) заменой g(t) на [...]; так как все произ- 
водные этой функции обращаются в нуль при #=8 и равны 
соответствующим производным функции (8 —?)*"* y(t) при 
==, то мы получаем выражение (12) для A,(x). Точно так 
же второй интеграл в правой части равенства (13) получается 
из интеграла (6) заменой y(t) на {...}; все производные 
{...} обращаются в нуль при f==a% и равны соответствую- 
щим производным (¢— a)" od) при ¢==8; поэтому из ра- 
венства (7) получаем выражение (12) для By(x). Тем самым 
формула (11) доказана. Если \====1, то из равенства (2) 
вытекает, что О(х-^) можно заменить на о (х-^\. 

Все наши результаты остаются справедливыми, если во 


всех формулах заменить на — 2, а в равенстве (7) выбрать 
путь интегрирования от # до — {<0 так, чтобы он целиком 
лежал в квадранте — п}2 < агв (и — я) = 0. Таким образом, 


получается описание поведения наших интегралов при x —> — со 
и тем самым асимптотическое поведение интегралов Фурье 
с тригонометрическими ядрами cos ХЁ и sin xt. 


2.9. Метод стационарной фазы 


Pacem отрим теперь интеграл 


Де lee at, (1) 


где х — большое положительное переменное и A(t) — веще- 
ственная функция вещественного переменного {. Согласно 
Стоксу и Кельвину, основной вклад в асимптотическое пове- 
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дение этого интеграла дают окрестности концов интервала 
и точек, в которых функция A(t) стационарна, то есть 
h' (1) =0; при этом, вообще говоря, вклад в главный член 
асимптотического разложения точек, где й’(1) ==0 более суще- 
ствен, чем вклад концов интервала. 

Предположим, что функция g непрерывна, а функция fh 
имеет непрерывную производную второго порядка; пусть 
<.— единственная стационарная точка функции h, «<*<. В, 
й' (<) =0 и 1’ (<) >0. Поскольку мы предполагаем, что 
основной вклад в значение интеграла дает окрестность точки 
7, введем новую переменную интегрирования и с помощью 
подстановки # (1) — h(t) = и*. Мы получим 


te 
Se~ J ge ata 


te 


= ( ои8 exp {ix {h (<) |} du, 


где u,==[h(t—s)—h(t)]2,. и, == [k(t +2) —1(*)]*. Так 
как лишь окрестность точки и==0 имеет существенное зна- 
чение, можно заменить g(t) Ha g(t) и выражение Зи/й’ (t) 
его пределом при #—+т, то есть выражением [2/й” (з)|"з. Мы 
получаем 
Ny 
у, : , 
29а] g(t) ( exp[ixu? + ixh (®) | du. 


uy 


В силу тех же соображений можно заменить интервал интег- 
рирования интервалом (— oo, со). Это приводит к формуле 


f~ [gra] Е екрив (9-Е жН], хо, (2) 


которая по существу и является результатом Кельвина. Вклад. 
в интеграл точки стационарной фазы t более важен, чем 
вклад концевых точек, поскольку с помощью интегрирования 
по частям можно показать, что если h’ (я) 520, й' (В) 520, 
то вклад концевых точек есть О(х-*). 

Принцип стационарной фазы применялся во многих матема- 
тических и физических задачах, но его трудно сформулировать 


3 А. Эрдейн 
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точным образом. По-видимому, наиболее пригодная теорема 
была дана Ватсоном (1920). Пуанкаре исследовал примене- 
ние принципа стационарной фазы к интегралам, содержащим 
аналитические функции; связь его работы с методом пере- 
вала указана в книге Копсона (1946). Метод стационарной 
фазы был также исследован Бийлом (1937) и в значительно 
более общем виде ван дер Корпутом (1934, 1936). 

Мы применим проведенное в предыдущем пункте иссле- 
дование интегралов Фурье, чтобы доказать теорему, которую 
можно рассматривать как точную формулировку и в то же 
время обобщение равенства (2). Назовем стационарной точ- 
кой порядка т, т=1, 2, ..., точку т, в которой A’ (<) = 
=h" (^) =... =” (3) =0 и A(t) 520. В окрестности 
такой точки A’ (1) == (Е—*)" А, (8), где №, (<) 520. Можно вве- 
сти и стационарные точки дробного порядка. Точка * назы- 
вается стационарной точкой (дробного) порядка в, если 
в некоторой окрестности этой точки A’ (f) имеет либо вид 
|1— |1, (№, либо вид (Ё—*)|Ё—*| 1, (1), где A, (<) 520. 
Если функция #(f) имеет на рассматриваемом интервале лишь 
конечное число стационарных точек положительного порядка, 
то можно представить интеграл в виде суммы конечного 
числа интегралов, в каждом из которых функция A(t) моно- 
тонна; не теряя общности, можно считать эту функцию в03- 
растающей. Поэтому мы ограничимся рассмотрением интегра- 
лов вида (1), в которых A(t) строго возрастает при «< 2< В, 
а аи В являются либо обыкновенными точками (то есть 
стационарными точками нулевого порядка), либо стационар- 
ными точками (положительного порядка). 

Пусть 0<), и=<1, 8( — функция, имеющая непре- 
рывные производные до N-20 порядка включительно на 
интервале а=<1=В, й(В — дифференцируемая функция, 
причем 


h(t) = (¢—a) (8 — 2)" 4, (0, (3) 


где р, в =1, и В, (1) — положительная функция на интер- 
sane «== В, имеющая непрерывные производные до N-20 
порядка включительно. Тогда 


8 
ее ¢ — a)" (B—AY dt = Bix) А, (4) 
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(где 
A(x) ~ Ay(x) и B(x) ~ By(x) до М-го члена при x — со, (5) 


a Ay (x) By(x) даются указанными ниже формулами (17) и (20). 
Введем сокращенное обозначение 


Г О Pe". (6) 


Для доказательства теоремы вновь используем нейтрализа- 
mop y(t), который является бесконечно дифференцируемой 
функцией на интервале a < = 8, такой, что при некотором 4, 
0<1< (8 —2)/2, имеем у (В == 1 при я = «-Рциу(й == 
при В —1={=3. Положим 
8-9 . 
— Ах) = { v(t) g, (eae, (7) 
a 
В 
ва = фи а, (t) =" at, (8) 
aty 
Чтобы получить асимптотическое разложение интеграла 
A(x), введем новую переменную интегрирования и, положив 


и = (И — 1 (5), ue h(B— 4) — h(a). (9) 


Из равенства (3) имеем 
1 : 
uw = h(t)——h(2) =) h' (s)ds = 


т 
1 


= {а fy ga ayy, (Е — yl dy, 


где S==a-+(t—a)y. Последний интеграл является положи- 
тельной возрастающей функцией от &, имеющей непрерывные 
производные до М-- 1-го порядка включительно. Поэтому 
равенство (9) задает отображение интервала & = #= 8 — 
на интервал 0 = и=и,, имеющее непрерывные производные 
до №М-- 1-го порядка включительно; следовательно, обратное 
отображение имеет непрерывные производные тех же порядков. 
Положим теперь у, (и) ==%(1) и : 


hu) = g(a, (10} 


3* 
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где 2, (1) дается формулой (6), а А (и) имеет непрерывные 
производные до N-ro порядка включительно на интервале 
0=и=и,. Тогда интеграл 

ay 

A(x) = — eine (y, Шиш" exp (ки) du 

° 
можно проинтегрировать по частям N pas, выбрав в качестве 
дифференцируемой части подынтегрального выражения у, А, а 
в качестве интегрируемой — все остальные сомножители под- 
ынтегрального выражения. Положим 


© 
Boner (t) = SE ( (2 — ura exp (6627) des (11) 
и 


результат интегрирования по частям можно записать в виде 


А (Хх) = Ак(х) + Юм(х), 
где 
N-1 F 
м(® = У (— 1" Е” (0) 9. „-, (0) (12) 


n=0 


Ry (x)= (— al они. (13) 


В интеграле (11) путем интегрирования является луч 
arg (2 — и) =п/(26) в комплексной плоскости. Очевидно, что 


(— It! pn fati т] ау 
P-n-1 0) = М г( 7 ) exo[ 58 | в. (14) 


Чтобы оценить Ф_„_, (и) при w>> 0, заметим, что |2 | -* == и” -* 


и 
ix’ +-x|z—u! p=ipx | ( --12—щехру, dt. 
9 


Так как вещественная часть последнего выражения отрица’ 
тельна, то 
[ехр (ix2") |< ехр(—х[2— и) 
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и, следовательно, 
lens) Se | |e—atrexo(—sle—al \d|z—ul< 
1\ 
<i (ats )e aes tad "+1, (15) 


С другой стороны, применяя метод перевала к интегралу (11), 
мы получаем, что при больших значениях xy’ 


1 (И) ==" О [(хи?)-"-1]. (16) 
Подставляя (14) и (15) в (12) и (13), имеем 
А, (х) == 
и 0) p(n ni(n +2) 
= pres ) [ : | = (mA O 2х8, 
У ag © Nos 7 OO Dp x е (17) 


a 


‘ -_! y(N м (а aN ый 
| Вл) [< Wop! (>) A [35 duX | 


Это показывает, что AWA, до N-ro члена, если hk < 1. Если 
\=1 и о=1, то тот же самый результат вытекает 
из п. 2.8. Пусть теперь № =1и p> 1; выберем $ так, чтобы 
выполнялось неравенство 

4 
aX, 2) 
du 


1 
(W= Tf ats i 


e 


1 
du<gé 


Так как из (16) вытекает, что Q_ y(t) = O(x~4) равномерно 
по и, когда и>8, то при достаточно больших x 


uy 


хм» Sie _yle) \| oe 


Таким образом, и в этом случае К==о(х-^). Тем самым 
доказан результат, касающийся A (x). 

Точно так же доказывается утверждение о B(x). Введем 
в формуле (8) новое переменное интегрирования, положив 


v= h(8)— A(t). (18) 


a du <. 
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Пусть 


Ko) =8, (QoS, (19) 
где g,(f) дается формулой (6). В повторных интегралах от 
<*-'ехр (—ixv’), возникающих при интегрировании по частям, 
мы интегрируем вдоль луча arg (2 — т) = —пй(25). Точно 
так же, как и для функции A, мы получаем, что B~ By, до 
№-го члена, где 


Ви — 


Ee = oe г ia i exp [ = hs ae | x ее ® (90) 


n=0 


В заключение этого пункта применим общий результат 
к функции 


1 
= | exp (ixt*) dt. 
о 


Здесь А=ЕН-==Т, р—8, o=— 1, = Ви, в силу 
выражения (17), 


Ак(х) = 1/8 Г(1/3} ей 8 x — 1, 


Далее, о=1— В, #==(1 — 9)"8, 


at i / 
19) = Baal —49)-* 
и 
Г 2/3 = 
ада Peay x) nai gfx 
по 
а потому 


SF ()NT (4/3) ert х- — у EP (oxy -n- те, ВБ. 
n=0 


Это разложение было получено в п. 2.7 с помощью метода 
перевала. Заметим, что в п. 2.7 х могло принимать комплексные 
значения, в то время как в нашем случае х—+ со, пробегая 
положительные значения. 
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ГЛАВА Ш 
ОСОБЫЕ ТОЧКИ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ 


В этой главе мы дадим краткое введение в асимптотичес- 
кую теорию обыкновенных однородных линейных дифферен- 
циальных уравнений второго порядка. Аналогичные теории 
существуют для уравнений любого (конечного) порядка и для 
систем дифференциальных уравнений первого порядка. Отно- 
сительно этих более общих теорий см. Айнс (1939, стр. 227 
и далее, стр. 645 и далее, стр. 653 и далее, стр. 685 и да- 
лее), Камке (1951, стр. 36 и далее, стр. 104 и далее, стр. 
131 и далее, стр. 160 и далее), Вазов (1953), литературу, 
указанную в этих книгах, и литературу, приведенную в конце 
этой главы. Асимптотические разложения встречаются также 
в связи с нелинейными дифференциальными уравнениями и 
дифференциальными уравнениями в частных производных. 

Мы изучим асимптотическое поведение решений уравнения 


У-Ер фу’ +(x) у=0 


при х—+х,. Здесь x либо вещественное переменное, пробе- 
гающее некоторый интервал (для которого х, обычно является 
одним из концов}, либо комплексное переменное, изменяюще- 
еся в некоторой области (для которой x, обычно является 
граничной точкой). В этой главе мы будем считать, что 
X,==00; очевидно, что это He приводит к потере общ- 
ности. 

Предполагается, что читатель знаком с основными теоре- 
мами существования для линейных дифференциальных урав- 
нений как в вещественной, так и в комплексной области, 
а также, что ему известны главные свойства решений этих 
уравнений. 
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3.1. Классификация особых точек 


В этом пункте мы изучим дифференциальные уравнения 
Ур у’ --9(Фу==0, (1) 


где х- комплексное переменное, пробегающее область R, 
задаваемую неравенствами т < | x] < <, a p(x) и g(x) — одно- 
значные аналитические функции в А (которые могут иметь 
особенности при х== <). Дадим краткий обзор хорошо из- 
вестной классификации изолированных особых точек уравнения 
(1) (см., например, Уиттекер и Ватсон, 1933, глава Х). За- 
метим, что в нашем случае особой точкой является х == со. 

Пусть у, (x) и у, (¥)— два линейно независимых решения 
уравнения (1), образующих фундаментальную систему; если 
аналитически продолжить эти функции вдоль некоторой кри- 
вой в Л, которая начинается и кончается в точке хи обходит 
в отрицательном направлении точку х==0, то мы получим 
две новые функции, которые можно обозначить через у, (xe ?™), 
J=1,2. Эти функции не совпадают с функциями у, (х), но 
они также являются решениями уравнения (1), а-потому имеют 
место соотношения вида 


у (хе-"") = a, y, (x) {-азу, (х), 
уз (107) = ау, (%) + eas (%)- (2) 
В этих соотношениях 
A= а; а,, 
Qs, as, 


— невырожденная матрица с постоянными элементами. 
Если заменить y,, у, другой фундаментальной системой, 
то матрица А заменится матрицей В, причем простое вычис- 
ление показывает, что B= МАМ-*, где М — невырожденная 
матрица с постоянными элементами. Таким образом, все ма- 
трицы, получаемые описанным выше способом, эквивалентны. 
Поэтому они имеют одни и те же собственные значения и 
одну и ту же каноническую форму. Поскольку собственные 
значения и кановическая форма этих матриц не зависят от 
выбора фундаментальной системы, они характеризуют особен- 
ность в точке х==с® (если эта особенность существует). 
Предположим, что собственные значения матрицы А раз- 
личны и она имеет поэтому диагональную каноническую 
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lox} 


Если y,, у, — фундаментальная система, соответствующая 
канонической форме матрицы А (каноническая фундаменталь- 
ная система), то соотношения (2) принимают внд 


yj (xe™) = у), f=, 2. (3) 


Положим Ху = exp (210) и назовем 2, 7, локазателялш, 
соответствующими точке со; каждое из чисел р,, 2, опреде- 
лено с точностью до целого слагаемого. Из равенства (3) 
следует, что каноническая фундаментальная система имеет вид 


форму 


y= x74, (x), j=l, 2, (4) 


rae 0, и &,— однозначные аналитические функции от x BR, 
которые, быть может, имеют особенности при x == со. 

Если собственные значения равны, }, =}, =) = exp (2112), 
то каноническая фундаментальная система имеет вид 


у, (x)= x кл (x), 
у. (x)= cy, ix) Inx +x7*d, (x), ) 


rae), и, нмеют те же свойства, что и в равенстве (4), 
а с-— постоянное; с==0, если каноническая форма матрицы A 
диагональна, и eset в противном случае. 

Точка х==ою называется обыкновенной точкой уравне- 
ния (1), если все решения регулярны в окрестности этой 
точки, то есть могут быть представлены в виде сходящихся 
степенных рядов 

а a, 
a ope aie aa es 

X == со называется регулярной особой точкой уравнения (1), 
если эта точка не является обыкновенной, а функции b, и, 
имеют полюсы при х==оо; в этом случае с помощью соот- 
ветствующего выбора показателей 2, и р, функции Ф, и $, 
могут быть сделаны регулярными на бесконечности; наконец, 
Х==со называется иррегулярной особой точкой уравне- 
ния (1), если по крайней мере одна из двух фуикций %,, $, 
имеет существенно особую точку на бесконечности. 

Можно показать (см., например, Уиттекер и Ватсон, 1933, 
стр. 278), что достаточным условием для того, чтобы 
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х = 00 было обыкновенной точкой, является 
P(x)== 2x7" Ох"), g(x) =0(“*, хо, (6) 


а достаточным условием для того, чтобы х==со было 
регулярной особой точкой, является 


P(x) =O(x"'), g(x) == O(x7*), хо. (7) 


В случае иррегулярной особенности ри 4 могут иметь суще- 
ственно особые точки при х==©; если ри 9 имеют лишь 
полюсы при х == 00, говорят об иррегулярных особых точках 
конечного порядка; наименьшее целое число А, для которого 


P(x) == 0(5^-'), g(x) == O(x**~*), хо, (8) 


называется рангом иррегулярной особой точки. Иногда регу- 
лярные особые точки рассматривают как особые точки нуле- 
вого ранга, 


3.2. Нормальные решения 


Если Х==<0 является обыкновенной точкой уравнения 
3.1(1), то у можно разложить в ряд по степеням x~'. Коэф- 
фициенты этого ряда можно вычислить с помощью рекуррент- 
ных соотношений, причем сам ряд сходится BR. Если x == 00 
является регулярной особой точкой уравнения 3.1(1), можно 
положить 


.7 


©. 
y= Dex", ¢, 0. 
n=0 


Мы получаем при этом квадратное уравнение для 0, рекур- 
рентные соотношения для с, (подобные приведенным ниже 
соотношениям (7)) и ряд для у, который сходится в КЮ. 
В обоих случаях коэффициенты могут быть легко вычислены, 
а сходящийся ряд удобно применять для вычисления значе- 
ний решения при больших х. 

Ситуация меняется коренным образом, если х==ео 
является иррегулярной особой точкой. Так как в этом слу- 
чае функции %,, Ф, имеют существенно особые точки при 
y= 00, TO следует положить 
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Мы получаем для с, бесконечную систему линейных уравне- 
ний, которую нельзя решить по рекуррентным формулам, а 
для р получается трансцендентное уравнение, содержащее 
бесконечный определитель (определитель системы). В этом 
случае коэффициенты трудно вычислить, а полученный ряд 
не будет быстро сходиться при больших значениях x, 

Томе открыл, что в случае, когда иррегулярная особая 
точка имеет конечный ранг, существуют формальные реше- 
ния, для которых описанные выше неудобства не имеют 
места; входящие в эти решения коэффициенты могут быть 
вычислены по рекуррентным формулам, а ряды достаточно 
удобны для вычислений при больших значениях х. Решения 
Томе имеют вид 


y = exp [P(x)] У Berle: ВО, 


п=0 


где Р(х) — многочлен; их называют нормальными реше- 
ниями. 

Мы покажем сейчас, как строятся нормальные решения 
в случае иррегулярной особой точки первого ранга. Заме- 
тим сначала, что если положить в уравнении 3.1(1) 


у==2ехр (—1, } рах), 


то для 2 получится дифференциальное уравнение вида (1), 
в которое не входит 2’. Таким образом, достаточно рассмот- 
реть дифференциальное уравнение 


У--9 (<) y=0, (1) 

в котором коэффициент g(x) задается сходящимся в А рядом. 
ao 

g(x) = У 4х". (2) 
n=0 


Мы будем искать формальное решение этого уравнения, 
имеющее вид 


2) 
у=е* У, Cx ee, ЕО, (3) 
n=0 


где хи р— постоянные. Следует отметить, что р в равен- 
стве (3), вообще говоря, не является показателем, принад- 
лежащим (в смысле и. 3.1) иррегулярной особой точке х == со. 
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При действиях над формальными рядами предполагаются, 
что 9-и„==0, ¢_,=0, m==1, 2,..., поэтому суммирова- 
ние производится от — со до -|- oo, за исключением случаев, 
когда это оговорено особо. 

Подставляя разложения (2) и (3) в уравнение (1), получаем 


o* Se,x7?-" — 20 J) (p +n) eno" -- 

EDGE EME M ege b Dagen® Doge t= 0. 

Сравнивая коэффициенты при одинаковых степенях, приходим 

к соотношениям 

we, — 20 (p +-n— 1)e,_, + (9 + a— 2)(p +-2—1)e,_,4 
+ У ae v=0, (4) 


yD 


справедливым при всех целых 7. Первое нетривиальное соот- 
ношение имеет место, если n=. Так как с, 50, получаем 


wo’ -++-9,=0. (5) 

Если в соотношении (4) 2==1 и © удовлетворяет равен- 
ству (5), то 

— 2up 4-9, =0. (6) 


Эти два уравнения определяют ® и р. Из равенств (4) можно 
получить также рекуррентные соотношения для коэффициен- 
тов. Заменяя в (4) п на п--1 и пользуясь равенствами (5) 
и (6), получаем 
2one, = (pm) (p--a— 1) e,_, + 
n-+-1 
ae > 4.6 +:-» П=12,... (Т) 
2 
Таким образом, нормальные решения существуют либо 
когда 4,50, либо когда 9,=4,==0. В первом случае 
равенство (5) определяет w, равенство (6) определяет р, а 
равенства (7) определяют коэффициенты (поскольку решения 
уравнения определены с точностью до постоянного множи- 
теля, можно считать, что с, =1). Кроме того, ®, р, С, ..., бт 
полностью определяются значениями Gy Qi, ... , Amar И 
обратно. В этом случае существуют два нормальных реше- 
ния, соответствующих двум возможным значениям ®. Во вто- 
ром случае, когда х == со является регулярной особой точкой, 
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из равенства (5) следует, что ®==0, равенство (6) обра- 
щается в тождество, равенство (7) при п==1 определяет р 
как один из корней квадратного уравнения p(p +-1)-+-9,==0, 
а равенства (7) при n==2, 3, ... определяют коэффициенты. 

Если 9, =0 и 4,520, то уравнения (5) и (6) несовместны, 
а потому нормального решения не существует. В этом слу- 
чае можно найти так называемое поднормальное решение. 
Для этого преобразуем заменой переменных 


=, 1(8) == "у (x) 


уравнение (1) в уравнение 
+ 3 
a+ [ее == = 1—0. (8) 


Если 9, ==0 и 9,520, то уравнение (8) имеет иррегулярную 
особенность первого ранга при х == со, а потому обладает 
нормальными решениями. В соответствии с этим при 9, =0, 
9,7£0 поднормальным решением уравнения (1) называется 
решение вида 


es} 
y = exp (x's) У, ogee hen, (9) 


n=0 


Относительно построения нормальных и поднормальных 
решений в случае особых точек более высокого ранга, а 
также для дифференциальных уравнений высших порядков 
и систем дифференциальных уравнений, см. Айнс (1939, 
стр. 653 и далее). 

Нормальные и поднормальные решения являются формаль- 
ными решениями; это означает, что если их подставить 
в дифференциальное уравнение и оперировать с ними, как 
со сходящимися рядами, то дифференциальное уравнение будет 
удовлетворяться. Однако бесконечные ряды, входящие в фор- 
мальные решения, вообще говоря, расходятся. Тем не менее, 
эти решения далеко не бесполезны, так как они дают асимп- 
тотические разложения решений уравнения (1). Связь фор- 
мальных решений с асимнтотическими разложениями решений 
была изучена многими авторами, начиная с Пуанкаре; неко- 
торые из этих работ указаны в конце главы. Следует отме- 
тить, что наиболее общие результаты получены Штернбер- 
rom (1920) для дифференциальных уравнений л-го порядка 
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любого (конечного) ранга и Тржитцинским (1933) для систем 
дифференциальных уравнений первого порядка. 

Вообще говоря, можно использовать два метода для 
доказательства того, что формальные решения являются 
асимптотическими разложениями решений дифференциального 
уравнения. Один из них восходит к Пуанкаре и был развит 
Горном в многочисленных работах, часть из которых указана 
в конце этой главы. Этот метод состоит в том, что для 
решений находят интегральные представления типа Лап- 
ласа, к которым применяют асимптотические разложения, опи- 
санные вп. 2.2. Второй метод, развитый Г. J. Биркгофом и его 
учениками, состоит в использовавии главных членов или ча- 
стичных сумм формальных решений для построения диффе- 
ренциального уравнения, в известном смысле близкого при 
болыших значениях х к заданному уравнению, и дальнейшем 
сравнении этих уравнений. В обоих методах важную роль 
играют сингулярные интегральные или интегро-дифференци- 
альные уравнения типа Вольтерра. 

Мы применим один из вариантов второго метода для 
исследования дифференциального уравнения (1) при 9,50. 
Приволимое ниже доказательство основывается главным обра- 
зом на работах Гогейзеля (1924) и Трикоми (1953, стр. 47—50). 
Так как аналитичность функций 9 и у не используется при 
построении формальных решений, исследование может быть 
проведено как в случае вещественного, так и в случае ком- 
плексного аргумента. 


3.3. Интегральное уравнение и его решение 


Рассмотрим сначала случай вещественного переменного, 
отложив краткое изложение случая комплексного переменного 
до п. 3.5. Рассмотрим уравнение 


y’ +(x) у=0, xpa>0, (1) 
где функция g(x) непрерывна при x >a 
$ a 
9%) > 9,87", x00, 9,0. (2) 
n=0 


Мы получаем в этом случае два формальных решения 


na 
em Sere", 6,560, (3) 
n=O 
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где @, р, с, ¢,, ... удовлетворяют соотношениям 3.25), 
(6), (7). Покажем, что эти формальные решения являются 
асимптотическими разложениями некоторых решений урав- 
нения (1). 

Пусть =, -|- №,; p—=p,-+-ip,; выберем ® = (—g,)"2 
так, чтобы функция е”^х-? была ограничена при х—+ со. 
Если 4, не является положительным вещественным числом, 
то этого можно достичь, выбирая значение квадратного 
корня, для которого ®, < 0; если 9, >0 и 10,520, то 
можно выбрать такое значение квадратного корня, что 9, = 
= Ве [4,/(2®)] >> 0; наконец, если 9, >0 и 9, — веществен- 
но, можно взять любое из значений квадратного корня. Таким 
образом, во всех случаях либо ®, > 0, либо ®, =0 и р, >0. 
В дальнейшем мы будем считать эти условия выполненными. 

Преобразуем уравнение (1) с помощью подстановки 


уе“). (4) 


Очевидно, что функция 2 удовлетворяет дифференциальному 
уравнению 


2-2 («—=) г-- 
+ [e ae + eet 1) Е 9 #=0. (5) 


В этом уравнении значения © и 2 удовлетворяют равенствам 
3.2(5), (6). Положим 


о) — 9, — 9.8") 4-89 ++ 1) == F (x). (6) 
Из соотношения (2) вытекает, что функция F(x) ограничена, 
|F(x)|<A, ха. (7) 


Перепишем уравнение (5) в следующем виде: 


а ~2p dz 5 - 
GE (etx "= ) Не х-22-2 F (x) 2 (x) =0. 
Интегрируя это равенство, получаем 

x 


d ваз ( tthe ee 
ete? ae a (О dt =e, x7 x", 
b 


где с; и b>>a— произвольные постоянные. Интегрируя еще 
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раз, приходим к уравнению 
x 
ик, t)F (t)2 (t)t-* dt =e, +e, \ e-*#at, (8) 
а 
где 


K(x, )=— ( exp [20 (f—s)] (+ ) * as. (9) 


м 


Уравнение (8) является интегральным уравнением типа 
Вольтерра. Любое решение уравнения (5) удовлетворяет при 
некоторых 6, с, и с, уравнению (8); обратно, любое решение 
уравнения (8) при любых 6,с,,с,, имеющее непрерывную про- 
изводную второго порядка, удовлетворяет уравнению (5). При 
р < сю существование таких решений уравнения (8) выте- 
кает из общей теории интегральных уравнений. Если же 
р = со, то (8) является сингулярным интегральным уравне- 
нием, а потому в этом случае надо доказать существование 
и дифференцируемость решений. 

Для того чтобы доказать, что уравнение (5) имеет реше- 
ние, которое может быть асимптотически представлено в виде 
ус„х“”, положим в уравнении (8) b==00, с, =1, с, =0. 
Уравнение (8) принимает при этом вид 


z(p= 14) K(x, YF eye tat. (10) 


Это интегральное уравнение можно решить методом последо- 
вательных приближений. 
Положим для любой функции &(х) 


Te (x)=) K(x, РОЗ tat (11) 


и определим далее функции 2,(х), 2,,,(x), 2(х) равен- 
ствами 

2, (х)=1, 2,., (х)== Tz, (x), п==0, 1, 2,..., (12) 

2 (х) = У, 2, (x). (13) 


п=0 
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Докажем, что функция 2(х) существует, удовлетворяет 
уравнению (10), дифференцируема и удовлетворяет уравне- 
нию (5). Доказательство разбивается на несколько шагов. 

Ядро K(x, t) ограничено при {2х7=х,, если x, >a 
и х, достаточно велико. 

Доказательство. Так как либо , < 0, либо ®, =0 
и р, 0, то при достаточно больших значениях $ имеем 


о 


а 53 2p, 
a п (е PSS?) = — Qo, SS >0. 
Следовательно, 
en PS 525, 


является возрастающей функцией от $. Положим 
tos) ($ \? 1-я [5 т 
ere (fs) (+) — et @-5) т [2 (5, ty + iv, (5, A], 
rae 


9, (s, )-+-is, (s, ty az ett (+). 


и применим к интегралу (9) вторую теорему о среднем зна- 
чении. Мы получим 
{ 


Ke, j= ( с. 4-9 (<)* lg, (s, f) +i, (5, 0145 = 


— et toy (4)" [ 


/ 


Я t 


+f ets, 045-29, (5, 945, 


5 % 


где х< 1= Интегралы в правой части этого равенства 
являются ограниченными функциями от x и Ё. Кроме того, 
если Ех 2х, и х, достаточно велико, то 


— x \? 
ее (2) ect, 
/ 


Поэтому 
[К(х. |< В, tax>x,, (14) 


при некоторых x, и В. 


3.3] ИНТЕГРАЛЬНОЕ УРАВНЕНИЕ И ЕГО РЕШЕНИЕ 83 


Если |5 (|< СЁ” при Ё2х,„, где ^ > -—1, To 


ABC nat 
ire + 


(15) 


| TE(x)|< 


В самом деле, из равенств (11), (7) и (14) вытекает 


[Te(x)| =] $ Kee NFS Ht tat < АВС "ЧИ, 


x 


откуда и следует неравенство (15). 

Для функций, определенных равенством (12), справедлива 
оценка 
(АВ 


т eh ож, (16) 


0 


|2. (x) [= 


Доказательство проводится по индукции. Оценка (16) сира- 
ведлива при и==0; если она справедлива при некотором м, 
то из определения 2„., и оценки (15) вытекает ее справед- 
ливость для п-|-1. 

Ряды (13) равномерно сходятся при x >> х, и функция 2 (x) 
удовлетворяет уравнению (10). Кроме того, 2(х) имеет не- 
прерывную производную второго порядка и удовлетворяет 
уравнению (5). Равномерная сходимость вытекает из оцен- 
ки (16). Если мы подставим = а, в интеграл, стоящий 
в правой части равенства (10). то, в силу равномерной схо- 
димости, можно выполнить почленное интегрирование. В силу 
соотношений (11), (12), (13), мы получаем, что функция z (x) 
удовлетворяет уравнению (10). Далее, интеграл в равен- 
стве (10) является дифференцируемой функцией oT x, а по- 
тому и функция z(x) дифференцируема. Так как 


io. OK (x, 1 jew t-x) [х \?? 
K(x, x)=0, Ox =e (+) ’ 
TO из равенства (10) вытекает 


[= 
z'(x) = fert-9 ("Feige tat, (17) 


х 
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Этот интеграл также является дифференцируемой функцией 
от x. Непосредственная подстановка показывает, что функ- 
ция 2(х) удовлетворяет уравнению (5). Если функция 2 за- 
дана формулой (13), то функция 


уе" (д) (18) 


удовлетворяет уравнению (1). В случае, когда 9, >0 и 4, 
вещественно, можно взять оба значения (— 9,)'' для ® и та- 
ким образом получить два линейно независимых решения 
вида (18). Во всех других случаях существует лишь одно 
решение этого вида, а второе решение записывается в виде 


га 


у, ху, (x) |, В, (19) 
b 


где р-—- любое число, достаточно болышое для того, чтобы 
при x 2b выполнялось неравенство у, (х) 5=0. Так как при 
х— с имеем 2(х)=1--О(х`'), такое $ заведомо суще- 
ствует. 

Таким образом, в любом случае на интервале x >= x, су- 
ществуют два линейно независимых решения уравнения (1). 
Если OX, Я эти решения можно распространить на интер- 
вал ха. Нам осталось показать, что полученные в п. 3.2 
формальные решения совпадают с асимптотическими разложе- 
ниями полученных в этом пункте решений. 


3.4. Асимптотические разложения решений 


Заметим сначала, что 


nm oO 

| e-tt-"at~en* У (—1)”")„х7’7”, когда x—> 00, (1) 

х m=0 ы 
г eS 
лее" ру ()„х-”-”, когда x —+ о, (2) 
5 то 


где 


(У, ==1, (У, =... и 1. (3) 
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Оба эти результата вытекают из формулы 2.1(5), если поло- 
жить в ней g==t~", й_„==(- 1)" ехр(-2 В. В частности, 


со 
ее О(е-*х-”, x + 00, (4) 
$ 


x 


el{- ‘df= O(e*x-"), хо. (5) 
b 


Докажем теперь по индукции, что функции, определен- 
ные формулой 3.312), могут быть разложены в acumnmo- 
тические степенные ряды вида 


зо 
2„(х) — 2, ий х— 00. (6) 


Это утверждение тривиально при #==0. Если оно верно для 
некоторого 7, то 


eS 
F(t)z,()~ 2 a,t-*, t—+ oo. 
Ho 
ы 2 
та = J em в-*( 5) ” F (t)z, (t) total. 


x 


Если фи р— вещественные числа, то, в силу последней 
теоремы п. 1.4, в интеграл можно подставить асимптотиче- 
ское разложение функции Р(1)2„(Ё) и проинтегрировать по- 
членно; при этом получаем 

Do ao 


, а 
У a § ewt-0 (+) НА x—+00. (7) 


kon x 


Если @ или р — комплексные числа, TO формула (7) может 
быть доказана путем подстановки в интеграл асимитотиче- 
ского разложения функции F(t) z, (f), состоящего из конеч- 
ного числа членов и остаточного члена, с последующей 
оценкой остаточного члена в формуле (7) по соотношению (4). 
Во всяком случае, каждый интеграл в равенстве (7) может 
быть разложен в асимптотический степенной ряд, который 
получается из формулы (1) и начинается с члена, содержа- 
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wero х-^-?. В силу третьей теоремы п. 1.4 можно подста- 


вить эти разложения в равенство (7). В самом деле, в (п) =, 
равномерность асимптотических разложений тривиальна, а 
ряды 1.4(5) состоят из конечного числа членов, так что 
вопрос об их сходимости не возникает. После указанной под- 
становки мы получаем путем перестановки членов, что 


Zaps ( (x) — = вых (8) 


Почлевное интегрирование разложения (8) приводит к иско- 
мому асимитотическому разложению функции 2,„.,. Наше 
утверждение доказано. 

В силу второй теоремы п. 1.4, можно подставить разло- 
жения (6) в формулу 3.3 (13). В самом деле, 2„==0О(х“”), 
равномерность разложений тривиальна, а ряды 1.4(2) конеч- 
ны, так что вопрос об их сходимости не возникает. Таким 
образом, функция z(x) может быть разложена в acumnmo- 
тический степенной ряд вида 


es] 
zix)~ Dc,x", x—> oo. (9) 


п=0 


\ 

Нам осталось доказать, что входящие вэто разложение 
коэффициенты с, удовлетворяют соотношениям 3.2 (7). 
Для этого заметим, что из ‚равенства 3.3 (17) и соответст- 
вующего соотношения для 2” вытекает возможность ‚разло- 
жения в асимптотические- степенные ряды функций 2” uz”, 
В силу результатов п. 1.6, отсюда следует, что разложение 
(9) может быть дважды почленно продифференцировано. По- 
лучающиеся асимитотические ряды должны формально удов- 
летворять уравнению 3.3 (5), а это и приводит к соогвошению 
3.2(7) для коэффициентов. Кроме того, с, =1. 

Мы доказали, таким образом, что функция у, (x), зада- 
ваемая формулой 3.3 (18), может быть асимптотически пред- 
ставлена одним из формальных решений; нам осталось пока- 
зать, что функция у,(х), задаваемая формулой 3.3 (19), 
представима с помощью другого формального решения. Для 
этого положим 


~ ex! 2, (Х). ^ (10) 


у: (®) = 
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Из равенств 3.3 (18) и (19) вытекает 


x 
2, (x) =2(x) | ero e-0( 2) (0-46 (11) 
b 
В силу выбора >, функция z(f) ограничена вне некоторой 
окрестности нуля; кроме того, как было показано выше, эта 
функция обладает асимитотическим разложением (9) при 
t— oo, где с, ==1. Поэтому при некоторых a, 
м-1 


(и = У, ae"tow-’, ть, 
n=O 
и 
2) N-1 x р is 
22 (X)__ 20 (x — t) Apes 
rar Dh OG + 


+) фене) О (Е- №) dt. 


Интегралы, стоящие под знаком суммы, можно асимитотически 
разложить по формуле (2), а последний интеграл, в силу 
оценки (5), есть О(х-^). Поскольку М является произволь- 
ным натуральным числом, 2, (х)/2(х) обладает разложением 
в асимитотический степенной ряд. Но тогда и функция 2, (x) 
может быть разложена в асимптотический степенной ряд 


2, (х) — ба х— <. (12) 


Точно так же, как и для функции 2, можно доказать, что 
входящие в это разложение коэффициенты удовлетворяют 
рекуррентным соотношениям, которые отличаются от соот- 
ношений 3.2(7) лишь тем, что ® и 2 заменены на —юи 
— 2. Таким образом, и решение 3.3 (19) асимптотически пред- 
ставимо одним из формальных решений. 

Если ® и о являются чисто мнимыми числами, то суще- 
ствуют два фундаментальных решения, представимых в виде 
3.3 (18); они определены однозначно с точностью до посто- 
янного множителя; оба эти решения ограничены и ни одно 
из них не стремится к нулю, когда х—+ oo. Во всех других 
случаях одно из фундаментальных решений, а именно 3.3 (18), 
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стремится к нулю, когда х-—> 00; это решение определено 
однозначно с точностью до постоянного множителя; второе 
решение 3.3 (19) не является ограниченным, когда х-—+ oo, 
Поскольку это второе решение зависит от 6, оно не является 
однозначно определенным. В самом деле, 


IS) FIIs (XH) ~ II, (4H), X00, YF. 


3.5. Комплексное переменное. Явление Стокса 


Результаты предыдущих пунктов могут быть распростра- 
нены на случай комплексного х, изменяющегося в области $, 
задаваемой неравенствами 


|| а a<argx<B. (1) 


Мы будем предполагать, что g(x) является аналитической 
функцией в 5, имеющей равномерно по arg x асимптотиче- 
ское разложение 3.3 (2), когда х—+ со в $. Мы имеем тогда 
два формальных решения 3.3 (3), где ® удовлетворяет урав- 
нению 3.2 (5). 

Линия Re [х (— 9,)**] =0 называется критической лини- 
ей или линией Стокса. Если х—* со вдоль одного из лучей 
критической линии, экспоненциальный множитель в обоих 
формальных решениях остается ограниченным, равно как и 
обратная ему величина. Если х—+ со вдоль любого луча, то 
главный член одного H3 формальных решений экспоненциально 
возрастает. 

Предположим сначала, что критическая линия не пере- 
секает область 5. Очевидно, в этом случае В —я< т, а 
потому существует решение ® уравнения 3.2 (5), такое, что 
Rewx < 0 при всех x в $. Если x изменяется вдоль любого 
луча argx==const в $, то результаты пп. 3.3 и 3.4 сохра- 
няют силу. Доказательство этих результатов для области $ 
протекает следующим образом. В интегральном уравнении 
3.3 (10) следует вести интегрирование вдоль луча arg х — ага 1. 
Тогда ограниченность ядра будет иметь место для любого x 
в 5 равномерно по x, если © < argx < В. Интегральное урав- 
ненне решается так же, как и выше, причем каждая функ- 
ция 2„(х) может быть выбрана Tak, чтобы она была анали- 
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тической в 5; при этом функция = (х) также аналитичнав $ 
как равномерный предел аналитических функций. В равенстве 
3.3 (8) 6 выбирается так, чтобы у, (х) 50, когда |x|), 
a<argx< 8. Отсюда вытекает существование двух реше- 
ний у, и у, в 5, асимптотические разложения которых отли- 
чаются лишь постоянными множителями от формальных реше- 
ний 3.3 (3). Асимптотические разложения имеют место, 
равномерно по argx, когда х—+ со, 2<argx<f. Любое 
решение дифференциального уравнения является линейной 
комбинацией решений у, и у,; асимитотические разложения 
этих решений получаются из асимптотических разложений 
У, и У». 

Рассмотрим теперь случай, когда критическая линия пере- 
секает область 5. Область $ распадается тогда на конечное 
число областей S,, R==1, 2, ..., К, отделенных друг от 
друга лучами критической линии. В каждой из областей 5, 
можно выбрать такое значение ®, для ®, что Rew,x <0 для 
всех x в 5,; таким образом, в каждом 5, существует фун- 
даментальная система решений V,,, Y,,, асимптотически пред- 
ставимых формальными рещениями. Относительно полного 
исследования этих фундаментальных систем см. Гогейзель 
(1924). Оказывается, что фундаментальная система решений, 
соответствующая лучу критической линии, может быть также 
сделана фундаментальной системой в каждой из двух обла- 
стей, разделенных этим лучом. Каждое из решений этой 
фундаментальной системы доминирует (экспоненциально воз- 
растает) в одной из двух областей и рецессивно (экспонен- 
циально убывает) в другой области. 

Рассмотрим решение у(х) дифференциального уравне- 
ния в 5. В каждой из областей $, у является линейной ком- 
бинацией двух фундаментальных решений, соответствующих 
этой области. Поэтому в каждой из областей у может быть 
асимптотически представлено как линейная комбинация двух 
формальных решений, однако коэффициенты этих линейных 
комбинаций могут меняться от области к области, Это обстоя- 
тельство было открыто Стоксом и называется явлением 
Стокса. Области $5, иногда называют областями (секто- 
рами) Стокса, а критические лучи — лучами Стокса. 

Относительно определения коэффициентов, входящих в вы- 
ражение у (x), как линейной комбинации формальных решений, 
см. Тьюриттин (1950). 
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3.6. Бесселевы функции нулевого порядка 


Проиллюстрируем результаты последних пунктов на при- 
мере дифференциального уравнения 


их" 2=0, (1) 
которому удовлетворяют функции Бесселя нулевого порядка. 


Замена переменной 
#=х-'у (2) 


преобразует уравнение (1) к стандартной форме 
” 1 с 
Уч о 


Это уравнение имеет вид 3.3(1), причем в неравенствах 3.5(1) 
мы можем положить а=0, aan В выбрать произвольно. 

Формальное решение этого уравнения получается так, как 
описано в п. 3.2; уравнения (5), (6), (7) этого пункта при- 
нимают в нашем случае вид 


w+1=0, p==0, 2wnc,=(n—1/2)*c,_,. 


Если ввести сокращенное обозначение a 


n 
То her (Pin + 1/2)? 
on I 2 2" nln (4) 
и выбрать соответствующим образом с 
шения уравнения (1) запишутся так: 


9» TO формальные ре- 


S, (x) = (2/m)"ne*— 4D) а, (— Nx" 'h, (5) 
n=0 
< 
$, (x) = (2/ny ee - А У отт", (6) 
n=0 


Между формальными рядами (5) и (6) имеют место тождест- 
венные соотношения 
5, (xe™) ==—S,(x), 5, (xe~™) ==5, (x), 
$, (xe™) = S, (x), $, (xe7™) == — $, (x). 


Так как ® ==-- 2 To критической линией является веце- 
ственная ось. В соответствии с теорией, изложенной в преды- 
дущем пункте, каждое решение уравнения (1) может быть 
асимптотически представлено в виде линейных комбинаций 


(7) 


3.6} БЕССЕЛЕВЫ ФУНКЦИИ НУЛЕВОГО ПОРЯДКА 91 


решений 5, и 5, в любом секторе, лежащем целиком в верх- 
ней (или в нижней) полуплоскости. При переходе через ве- 
щественную ось коэффициенты этих лннейных комбинаций 
могут изменяться. Мы увидим сейчас, что такое изменение 
действительно имеет место. 

Непосредственная подстановка показывает, что одним из 
решений уравнения (1) является бесселева функция первого 
рода нулевого порядка 


ти ат 

лы=У = (s+) , (8) 
m=0 

которая является четной целой функцией oT x. С помощью 

разложения показательной функции в степенной ряд и по- 

членного интегрирования легко убедиться, что для функции 

J,(x) справедливо интегральное представление Пуассона 


J, (x) == хи (1 — и)" du. (9) 


Предположим, что Rex > 0. Тогда e' экспоненциально 
убывает, когла ти—» <, а потому интеграл (9) можно 
представить в виде 


| =- 


В первом интеграле положим 


a==1-+it, фи", | 4-u= 2+ it,. 


thin РЕГ 


а BO втором положим 
и=-А-И, 1—u=2—it, Ти, 


Мы получаем при этом две функции, отличающиеся лишь по= 
стоянными множителями от функций 


wo 
НР == ава | ее) 4 
. (1 0) 
BE (x)= елена { en 8g" (2 — it) dt, 
$ 


92 ОСОБЫЕ ТОЧКИ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ [Гл. Ш 


Функции, определяемые равенствами (10), называются бессе- 
левыли функциями третьего рода (или функциями Ган- 
келя) нулевого порядка. Эти функции определяются равен- 
ствами (19) при Кех_>0. Однако область их определения 
можно распространить до —п< ая х< 2m в случае но 
и до —2тп<аех<п в случае Н®. Для этого надо 
повернуть путь интегрирования так, как это было описано 
вп. 2.2. 

Можно показать, что № и Н® также являются реше- 
ниями уравнения (1). Очевидно, что 


Jy) = HY (+E HP (09, (11) 


причем, как показывает более тщательное исследование, обе 
функции HY и HY имеют логарифмические особенности при 
х==0. Изучение этих особенностей дает возможность опре- 
делить функции Ганкеля для всех значений arg x. 

Интегральные представления функций Ганкеля являются 
интегралами Лапласа, а потому их асимптотические разложе- 
ния при больших значениях х могут быть получены с по- 
мощью последней теоремы п. 2.2. Мы получаем при этом, 
что 


НР (x) ~ S, (x) (12) 

равномерно по arg x, когда х — oo, —п-- = агех < жт— в, 
=> 0. 

HY? (x) ~ 5, (x) (13) 


равномерно по arg х, когда x —+ со, —2n-Le<argx<n—e, 


=> 0. 


Из равенства (11) вытекает 
QJ, (x) — $, (x) +S, (x) (14) 


равномерно по arg x, когда х—+ со, —п-|-= = а19 хп — в, 
=> 0. 

В последнем равенстве символ ~ применен в несколько 
необычном смысле, поскольку в правой части равенства мы 
имеем не одно асимптотическое разложение, а сумму двух 
таких разложений. Оправданием такого использования знака — 
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является то обстоятельство, что на вещественной оси оба 
разложения имеют один и тот же порядок, и сумма (14) 
может быть записана в виде одного разложения (которое 
уже не является асимптотическим степенным рядом), в то 
время как в верхней (соответственно в нижней) полуплоскости 
слагаемое S,(x) (соответственно $, (х)) рецессивно и может 
быть опущено. 

Мы получили, таким образом, асимптотическое разложе- 
ние функции J, (x), справедливое во всей плоскости, за исклю- 
чением узкого сектора, окружающего отрицательную полуось. 
Чтобы получить асимптотическое разложение, справедливое 
и в секторах, содержащих эту полуось, заметим, что из 
равенств (8) и (14) вытекает соотношение 


QJ, (x) = 2, (xe™) — $, (хе) +S, (хе) 


при x—+co nH —п-- == ас (хе) =п—е. Поэтому, в силу 
равенства (7), 


2J, (x) ~ $, (x) — 5, (x) (15) 

равномерно по arg х, когда х — ©, —2п +e <Sargxa—e. 
Аналогично, J,(x)==J,(xe~™) и 

2J, (x) ~ — $, (x) +S, (x) (16) 


равномерно по argx, когда х—+ co, е= ао х = 2т— в. 

Сравнение равенств (14), (15) и (16) показывает, что 
в нашем случае имеет место явление Стокса. Лучи, при- 
надлежащие узким секторам, являются лучами Стокса. На 
первый взгляд может показаться странным, что секторы, 
в которых справедливы различные асимитотические разложе- 
ния, перекрываются друг с другом, однако в этом нет ника- 
кого противоречия. Например, области, в которых справедливы 
разложения (14) и (15), имеют общую часть 


— пах = —е. 


В этой общей части слагаемое S,(x) рецессивно, а потому 
правые части равенств (14) и (15) асимптотически равны. 
Таким образом, коэффициенты при формальных рядах скачко- 
образно меняются в секторах, где эти ряды доминируются 
другими рядами, 


94 ЛИТЕРАТУРА 


ЛИТЕРАТУРА 
Айнс 3, Л,, Обыкновенные дифференциальные уравнения, Харь- 


ro. faa 

Камке ' Э.' Справочник по обыкновенным дифференциальным урав- 
нениям, Москва, 1951. 

Униттекер Е. Т. и Ватсон Г. Н., Курс современного анализа, 
Л. —М., 1933. 

Hoheisel О. К. H., J. Reine Angew. Math. 153, 228-—248, 1924. 

Horn Jakob, Jahresber. Deutsch. Math. Verein. 24, 309—329, 1915. 

Horn Jakob, Jahrester, Deutsch. Math. Verein. 25, 74—83, 
301—325, 1916. 

Horn Jakob, Math. 7. 3, 265—313, 1919. 

Horn Jakob, Math. 2. 8, 100—114, 1920. 

Horn Jakob, Math. Z. 21, 85—95, 1924. 

Poincaré Henri, Acta Math. 8, 295—344, 1886. 

Sternberg Wolfgang, Math. Ann. 81, 119—186, 1920. 

Tricomi Е. C., Equazioni differenziali, Second edition, Torino, 1953. 

Trjitzinsky W. J., Acta Math. 62, 167—226, 1933. 

Turrittin H.L., Trans. Amer. Math. Soc. 68, 304—329, 1950, 

Wasow W.R., Introduction to the asymptotic theory of ordinary 
linear differential equations. Working paper National Bureau of 
Standards, 1953. 


ГЛАВА IV 


ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ 
С БОЛЬШИМ HAPAMETPOM 


В этой главе будет кратко изложена асимптотическая 
теория обыкновенных однородных линейных дифференциаль- 
ных уравнений второго порядка, содержащих большой пара- 
метр. Относительно этой теории, а также ее обобщений на 
уравнения высших порядков и на системы дифференциальных 
уравнений первого порядка, см. Айнс (1989, стр. 365 и далее), 
Камке (1951, стр. 74 и далее, стр. 111 и далее, стр. 163 
и далее, стр. 209 и далее), Вазов (1953), литературу, ука- 
занную в этих работах, а также литературу, указанную 
в конце этой главы. 

Как и вп. 3.2, дифференциальные уравнения можно 
преобразовать к стандартному виду 

УЧ, Л у=0, 
где х— вещественное или комплексное переменное, а ) — ве- 
щественный или комплексный параметр. Мы изучим поведе- 
ние решений этого дифференциального уравнения при h—> dj, 
причем, не теряя общности, можно считать, что }, = со. 

Предполагается, что читатель знаком с основными теоре- 
мами о зависимости решений дифференциального уравнения 


от входящих в это уравнение параметров, см., например, 
И. Г. Петровский (1939, стр. 54). 


4.1. Проблема Лиувилля 


В ходе своих классических исследований по проблеме 
Штурма — Лиувилля МЛиувилль изучил поведение решений 
дифференциального уравнения 


Ур) +r (x)] у=0 (1) 
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при \—+ oo. Здесь x — вещественное переменное, а < х =, 
р(х) — положительная функция, имеющая непрерывную вто- 
рую производную, а г(х) — функция, непрерывная на интер- 
pane а==х=6. Метод Лиувилля может быть описан сле- 
дующим образом. 

Введем новые переменные $ и 1 с помощью подстановки 


= ро] ах, перу. (2) 


Ora подстановка преобразует интервал а == х = в интервал 
9 == 8, а дифференциальное уравнение (1) — в уравнение 


р (3) 
rae 
es Lp? 5 р"? r 
ив ть (4) 


Pp 
является непрерывной функцией or 6, 2<E<B (р’==ар/ах 
ит. д.). 
Рассуждения, аналогичные проведенным в п. 3.3, показы- 
вают, что решения уравнения (3) удовлетворяют интеграль- 
ному уравнению Вольтерра 


В 
(==, cos 38 с, чп 4-Е №-* | sink —0р(01(046 (5) 


т 


где а == в, а C,, 2, — произвольные постоянные. Функции 
1 (5) ис, соз Е -- с, sin k§, равно как и их производные, имеют 
при &==7 одинаковые значения. 

Применяя метод последовательных приближений, можно 
получить решение интегрального уравнения (5) в виде 


(6) 


где 
ny (Е, 4) ==, cos XE с, sin 24, 


Е 
Чень & МА" ( sink E—A oH ga Nat, m=O, т... 
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Если |2()|=< А, то легко доказать по индукции, что 


еее» "Е" | 
16 “al - И al g ed, Dy cus 


|» G ^) [= 


В случае, когла интервал (%, 8) конечен, отсюда вытекает, 
что ряд (6) равномерно сходится при & ==, №2), >0, 
а также, что он является асимптотическим разложением функ- 
ции 1(5, 4) при h—+ ce. 

Функции %, довольно трудно вычислить. Поэтому при 
больших значениях ). используют другие приближения, по- 
лучаемые из формальных решений (которые, вообще говоря, 
расходятся). Для этого существуют два метода. Применяют 
либо формальные разложения 


> в, (x) 47" cos. § 4- У 8, (x) 47" sin АЕ 


для у(х, }), либо формальные разложения 


Dn 
aE+ Ху," 

п=о0 
для пу(х, №). Во втором случае у является решением, не 
имеющим нулей. В обоих случаях приближения для у строятся 
с помощью частичных сумм формальных разложений; эти 
приближения сравниваются с у с помощью интегрального 
уравнения. Если частичные суммы состоят из одного члена, 
то оба метода сводятся к процессу Лиувилля. 


4.2. Формальные решения 


Вместо уравнения 4.1(1) мы будем рассматривать несколько 
более общее дифференциальное уравнение 


+x, Му=о. (1) 


Если g(x, Х) является формальным степенным рядом, распо- 
ложенным по степеням 47", коэффициенты которого зависят 
от X, TO два линейно независимых решения уравнения (1) 
также могут быть представлены в виде формальных степен- 
ных рядов по }-'. Если же формальное разложение функ- 
ции 9 по степеням \ содержит положительные степени }, то 
формальные разложения у имеют вид ряда Лорана. Тем 
не менее мы увидим, что в случае, когда q(x, №) является 


4 А. Эрдейи 
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функцией от 2, имеющей лишь полюс при A== со, можно 
построить формальные решения, аналогичные нормальным и 
поднормальным решениям из п. 3.2. 

Предположим, что в уравнении (1) g(x, h) имеет вид 


< ` 

Хи", (2) 

n=0 
где функции 9„(х) независимы от 2, а # — положительное 
целое число. Предположим также, что функция 9,(х) не 
обращается в нуль в интервале или односвязной области 
комплексной плоскости, где изменяется переменное х. 

В соответствии с двумя методами, указанными в конце 

п. 4.1, существуют два типа формальных решений уравнения 
(1). Первый из них имеет вид 


a : Rk-1 
D 2%) A" exp [ SH 8, (x) MY. (3) 

n=0 v=0 
Подставим разложение (3) в уравнение (1); при этом 
условимся считать, что g,==0, 9„==0 при n==—I, —2, 


—3,...,4 В, =0 при у==—1, —2,... и при у==А, #--1,... 
Тогда суммирования можно считать распространенными на 
все целые числа и потому 


мм ед" 
+2DH MDa AM Warn aa" Sa dn = 0. 


Приравнивая нулю коэффициенты при 4°47", получаем, что 
при всех целых значениях и выполняются равенства 


ии, Е 
+2 Ха, Pn etn = 0- (4) 


Первое нетривиальное условие соответствует п==0. Если 
положить в (4) п=0,1,..., &—1, то получим 


„НУВ, yO =O; Ay vey Rd 
ИЛИ 


Ви =0, (5) 
Оч У 0, и... 6) 
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Вследствие этих соотношений в первой сумме равенства (4) 
можно ограничиться суммированием по 2 А. Если положить 
в равенстве (4) л==А, то придем к равенству 

R-1 


2a) 8 +45 (By +e + У BR, =O (7) 


У=1 


а если заменить в равенстве (4) п на &--п, то получим 


Е: 


20, By tn (Bo at У Ви) 


v= 
п Ат 
С =P >» By Bet my) 
та УЕ: 
п 
--2 Ух, Bite, ,=0 п=1, 2,... (8) 


тт 


Мы видим, таким образом, что если а, и В, удовле- 
творяют соотношениям (5) — (8), то ряд (3) формально удо- 
влетворяет уравнению (1). При этом пустые суммы (то есть 
суммы, верхний предел суммирования которых меныше ниж- 
него предела) считаются равными нулю. Так как 49,520, 
то можно выбрать одну из ветвей функции |-— 4, (х)|"*. Тогда 
равенство (5) определяет 8, с точностью до произвольного 
постоянного. Кроме того, 8 £0 и, следовательно, равен- 
ства (6) рекуррентно определяют функции §,,...,8,_,, каж- 
дая из которых также определена с точностью до произволь- 
ного постоянного. Уравнение (7) позволяет найти a, с точ- 
ностью до постоянного множителя, после чего из уравне- 
ния (8) можно рекуррентно найти функции 9,, 2, ..., каждая 
из которых определена с точностью до слагаемого вида Са, 
где С, — постоянная величина. Соответственно двум ветвям 
(—4q,)'2 мы получаем два формальных решения, имеющих 
вид (3). 

Второй тип формальных решений имеет вид 


© rs 
exp | >) В, (x)a*"]. (9) 


n=0 
Подставляя это выражение в уравнение (1), имеем 
star, А-н ур’ ЗА-м? a te-n 
a Pak of (УВ, м") Уча =0, 
4* 
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Приравнивая нулю коэффициенты при }**-", получаем соот- 


ношения 
В --4%=0, (10) 
п—1 
28,8. dn-t Я М w= A=1,...,8—1, (Il) 
п =} 
23,8, 9, - i el ==0, n=k, k+1, suet, (IQ) 


Существуют два линейно независимых формальных решения 
этого тина. 

Связь между этими двумя типами формальных решений 
совершенно очевидна. Уравнения (10) и (11) тождественны 
с уравнениями (5) и (6), а Ух,\-" является формальным 
разложением 


exp( Saat) в 


nak 


На протяжении этого пункта мы предполагали, что 
q(x, k) как функция от h имеет при ).=с<®о полюс четного 
порядка. Если полюс имеет нечетный порядок, то решений 
вида (3) и (9) не существует, и вместо разложения по сте- 
пеням A надо разлагать по степеням 4/2, 


4.3. Асимптотические решения 


Мы покажем теперь, что при некоторых предположениях 
дифференциальное уравнение 4.2(1) обладает фундаменталь- 
ной системой решений, которые асимитотически иредставля- 
ются полученными в предыдущем пункте формальными ре- 
шениями. При этом несущественно, будем ли мы сравнивать 
решения уравнения 4.2(1) с выражениями 

№- й А-1 . 
Sa, (2) kM exp |S B(x) MY], 
n=0 У=0 
где a, и В, уловлетворяют соотношениям 4.2(5) — (8), или 


с выражениями 
+N-1 


exp "5 м-р 


и=о 


4.3] АСИМПТОТИЧЕСКИЕ РЕШЕНИЯ 101 


где В, удовлетворяют соотношениям 4.2(10)— (12). В самом 

деле, коэффициенты & и В можно выбрать так, чтобы OTHO- 

шение этих двух выражений имело вид 1--О(\-^). 
Зафиксируем положительное целое число N и положим 


a #-т у 
Y¥,(x)=exp[ > ВМ", 7=1, 2, (1) 
n=0 
где BY —=— 8), и при каждом / коэффициенты В,, удовле- 


творяют соотношениям 4.2(10) — (12). Эти коэффициенты 
полностью определяются функциями Gy, ...,G2e4+N—1 И не- 
которыми производными этих функций. Будем говорить, что 
функции J, имеют достаточное число производных, если все 
производные этих функций, нужные для определения функ- 
ций Вар п=0,...,2#--№М— 1, существуют и непрерывны. 
Мы считаем, что х изменяется в ограниченном замкнутом 
интервале /: a= х==р, a) изменяется в области 5: || 2=)., 
ф, == аго) == $,. Теорема, которую мы хотим доказать, форму- 
пируется следующим образом: 

Пусть для любого фиксированного \ в S q(x, i) яв- 
ляется непрерывной функцией от x на интервале I; пусть, 
далее, 


ЗАЕМ -1 
q(x, y= р Gn (x) \** "+ 00-^) (2) 


равномерно по x и argh, когда \—+ со в $. Здесь 4,(х)— 
функции, дифференцируемые достаточное число раз на 
интервале I, причем 


Re {1 [—g, (x)]'2} 40, (3) 


если \ принадлежит S и x принадлежит 1. Тогда диффе- 
ренциальное уравнение 


y+ a(x, d) у==0 (4) 

имеет такую фундаментальную систему решений у, (х), 
у» (х), что 

= И -- 00-м), 6) 


=, [14+ 00-%)] 


равномерно по x и argh, когда }—w в $, 
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Докажем эту теорему методом, аналогичным использо- 
ванному в п. 3.3. В силу соотношений (3) и 4.2(10) можно 
выбрать В: и By, так, чтобы для любого h в $ Ве [h*8,, (x)] 
была возрастающей функцией от x, а Re [h*8,, (x)] — убы- 
вающей функцией. Из равенства (1) следует тогда, что для 
любого достаточно большого } из области $, |У, (x)| являет- 
ся возрастающей, а |У,(х)| — убывающей функцией от x. 

Чтобы установить существование и асимптотические свой- 
ства решения у, (x), сделаем в уравнении (4) подстановку 

у, (x) = У, (x) 2 (х). (6) 
Мы получим 


т 
Zz +22 -+- F(x, №) 2==0. (7) 
1 
rie 
у’ 
F(x, == у, 4-9 = 
2ktN-1 2k-+N~-1 
a tO pe vem tg. (8) 
n=0 n=0 
В силу соотношений (2) и 4.2(10) — (12), имеем 
F(x, )=00-^ 


равномерно по x и argh, когда h-> co в 5. Уравнение (7) 
можно записать в виде 
а 


ия | + и 9 2=0. 


Дважды интегрируя и выбирая соответствующим образом про- 
извольные постоянные, получим : 


z(x)=1— K(x, 0 Fit, Nz Hat, (9) 


где 
х 


K(x, = У: "(9 45. 
t 


Так как функция |У, (х)| возрастает, то | У, (| == | ¥,(s)| и 
А, |= 0—9, == х=5. 
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Существование 2(х) вытекает теперь из общей теории интс- 
гральных уравнений Вольтерра либо может быть установлено 
< помощью метода последовательных приближений. Из ра- 
венств (8) и (9) вытекает, что когда } —› сов 5, то 2(х) = 
=—=1--О(-^) равномерно по x и argh. Кроме того, функ- 
ция 2(х) дифференцируема, причем 

х 

ИУ PW) F U,V) z() dt = O00, 

a 

Но тогда 


У, (<) = 
=, (x) [29 a г | = ¥,(x)[1-+ 00.-^)|, 


(x 


Тем самым равенство (5) доказано при /==1. Доказа- 
тельство при /==2 протекает точно так же, за исключе- 
нием того, что в качестве постоянного предела интегриро- 
вания в интегральном уравнении принимается не а, а 6. 


4.4. Приложения к бесселевым функциям 


Мы применим теперь методы, развитые в предыдущих 
двух пунктах для доказательства асимптотических формул *) 


J, 0 зв 8) ~ (27) th В) "2 exp (kthB — 7.3), Х-+ою, (1) 


Hi” (sec В) ("> tg в) 72° ехр[-- 1488 — 8 — п/4], (2) 
h— со. 


Равенство (1) имеет место при 8>>0 равномерно по В в ин- 
тервале 0 < В, <3 <8, < со. Равенство (2) справедливо при 
0<3< 1/2, равномерно по 3 в интервале 0 === 1/2 — =; 
в этом равенстве (а также ниже в равенстве (18)) верхний 


знак соответствует функции He, а нижний знак — функции 


Н®. Оба эти результата могут быть выведены из интеграль- 
ных представлений бесселевых функций с помощью метода 
перевала (Ватсон, 1949, пи. 8.4 и 8.41). 


*) Не следует смешивать использованный в этих формулах сим- 
вол } с используемой ниже функцией } (x). 
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Функции 
DCX), хЪу, 0х), cH? (hx), x? (dx) 
являются решениями дифференциального уравнения 
y +P — 0 — 1/4) х-] у=0. (3) 
Это уравнение имеет вид 4.3 (4), где 
k=1, 9, (5х) =1—х7*, 9,(х) ==(2х)", 


а остальные 4, (х) тождественно равны нулю. Точки х==0, со 
являются особыми точками уравнения (3), а х=1 является 
так называемой точкой перехода, то есть точкой, в которой 
при всех значениях \ нарушается условие 4.3 (3). На любом 
интервале а=х==6, не содержащем ни особых точек, ни 
точки перехода, теорема из п. 4.3 дает общий вид асимпто- 
тических решений, но она недостаточна для того, чтобы по- 
лучить выражения у, и у, через стандартные бесселевы 
функции. Для того чтобы выразить наше решение через бес- 
селевы функции, увеличим интервал так, чтобы один из его 
концов совпал с особой точкой уравнения (3). В этом случае 
уже нельзя воспользоваться теоремой п. 4.3; однако, приме- 
няя использованные там методы, мы получим требуемый ре- 
зультат. 

Рассмотрим сначала уравнение (3) на интервале 0 =х OCI, 
Из равенства 4.2 (5) вытекает 


В =k (x? — 1)" 
и, следовательно, 
ВЯ Пн 
= (Е — ха 
Из 4.2 (7) получаем 


i а= =). (4) 
20,8, +48) =0 
и потому 


а, (x) == [В' (x)] - V2 =o (x) == (1-х?) 7, (5) 


С помощью определенных нами функций & и В образуем 


функции 
У, =a (x)e'P, У, (x)= a (хе, (6) 


которые соответствуют главным членам формальных решений. 
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Интегральное уравнение для z==y,/¥, имеет вид 


(к) = Ь— КР 9 2(0 46 (7) 
о 
Непосредственное вычисление показывает, что 
ое 30 44x 
reaper (id) dette, © 


а потому функция F ограничена на нащем интервале. Кроме 
того, из равенства (5) следует, что х7*==8’ и потому 


‚ОГ 0 0 фм 
t 


ot — exp {2% [8 (9) —8 (>. (9) 


Ho B(x) является возрастающей функцией, В (8) — B(x) <0 
пр 0<Ё=х=<6< 1. Поэтому  показательная функция 
в формуле (9) ограничена, если Ке» 0; ограничена и 
функция a? (1). Мы получаем, таким образом, оценку 
FU Kw, |= ТЕ, (10) 

справедливую при Reh0, Ota x<b<1. Коэффициент 
С не зависит от }h, хи Ё. 

Из полученных результатов легко вытекает, что интег- 
ральное уравнение (7) можно решить с помощью метода по- 
следовательных приближений, причем решение имеет вид 


zx) ¥ z, (x), 


п=:9 


где 
2, (х) =1, 
x 


ти (= | K(x, NF 2, 2=0, 1... 


В самом деле, из оценки (10) легко получить по индукции 


ину (Gir) + 


106 ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С БОЛЬШИМ ПАРАМЕТРОМ [TJ IV 


Поэтому ряд, определяющий 2, сходится равномерно по хи 
\, если Х лежит вне некоторой окрестности нуля; функция 
2(х) удовлетворяет интегральному уравнению, имеет непре- 
рывную вторую производную, удовлетворяет соответствующему 
дифференциальному уравнению И такова, что 


2=1--0(+). 


Тем самым установлено существование такого решения у, 
уравнения (3), что 


у, [1-0 (4), и 
O<e<d<1, ВА. 


Нам осталось показать, что это решение отличается лишь 
постоянным множителем OT J) (AX). 
Так как у, является решением уравнения (3), то 


ху, (x) ==с, (h) Л 0х) +c, (h) Л 0х). (12) 


Зафиксируем теперь * и рассмотрим равенство (12) при x —+0. 
Хорошо известно (Ватсон, 1949, стр. 55 и 78), что 


Из соотношений (11), (6), (5) и (4) вытекает, что *) 
May, (X)~ oY, (x) ~ хр eh, 0, 


Поэтому, разделив обе части равенства (12) на (х/2)’ и пе- 
рейдя к пределу, при х—+0 получаем, пользуясь указанными 
выше соотношениями, что 


(= -ТО-НТ, с,0)=0. 


По формуле Стирлинга отсюда следует 


c, 0) == (2 [1 +0 (+)]. 


*) Заметим, что У, — приближенное решение, У, — бесселева 
функция второго рода. ‘5 
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Мы доказали, таким образом, равенство 


Аи 
пе у, (x) [1-+0(= Wee 
= (тах) - 1+ ¥, (x) [1 +o(+)] , (13) 


где Ох <! и Ве} =0. Если положить х==зес 3, 
подставить вместо Х выражение (6) (принимая во внимание 
равенства (4) и (5)) и считать } в последней формуле поло- 
жительным, мы получим равенство (1). 

Рассмотрим теперь уравнение (3) на интервале 1<а-х<« оо. 
В этом случае 


а (x) == x" (x? — 1) > == [-— 18" (x) -*, (14) 
B(x) =i \ (зах = i (x? — 1)": — Рагссоз x7! 


где arccos х — главное значение функции, обратной функции 
cos x, в частности, агссоз (х-') —* п/2, когда х—+ со. Pynk- 
ции сравнения имеют теперь вид (6), где a и В определяются 
формулами (14). Интегральное уравнение для 2==у,|У, имеет 
вид 

oO 


z= 14+) Ах, DPE Wz at. (15) 


Равенство (8) сохраняет силу и в этом случае; оно показы- 
вает, что F(t, 4)== O(t-*) при #26531 и всех ). Так как 
в нашем случае a~*==—8’, то вычисление интеграла в фор- 
муле (9) приводит к равенству 


2 (1 — exp {24 [80—82]. 


Но — 3 является возрастающей функцией от x, а Ё=х 
поэтому показательная функция в полученном равенстве огра- 
ничена, если Imh >20. Ограничена и функция a” (В), а потому 
оценка 


K (x, t) = 


C 
ile 


РВК, A < Т (16) 


справедлива при Imk>0, 1«<я=х=й коэффициент С 
в этой оценке не зависит OT h, x, ft 
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Положим теперь 


25 (Х) =1, 


Zar (=) K(x, NF 2, dt, п==0,1,..., 


2x) р %,, We). 


Из оценки (16) легко вытекает по индукции, что 
ГС owe 
<— (| 
|2, (9) | = п! (=) : 
а потому ряд, определяющий 2, равномерно сходится, если } 
лежит вне некоторой фиксированной окрестности нуля. Функ- 
ция 2(х) удовлетворяет интегральному уравнению, имеет 
непрерывную вторую производную, удовлетворяет соответ- 
ствующему дифференциальному уравнению и такова, что 


ад=1--0 (3). 


Отсюда вытекает существование такого решения у, уравне- 
ция (3), что 


У: (x) = 7, (%) [! =F 9(5+) |, (17) 
1<а=х«®, ш\ 20. 


Нам осталось показать, что это решение отличается от 


HE (ix) лишь постоянным множителем. 
Так как является решением уравнения (3), то 
1 р р: > 


xy, (x) =, () AM (hx) -- с, 0) Н® Ox). (18) 


Зафиксируем значение \ и рассмотрим поведение обеих ча- 
стей этого равенства при х—> oo. Хорошо известно (Ватсон, 
1949, п. 7.2), что при 0 < аге < п имеем 


5 2 \'s р d 
HO A, x) ~ (==) * ехр [= > z)]. х— 00, 
а из соотношений (17) и (14) вытекает, что 


у, (x) ~ У, (x) ~ exp [18 (x)] ~ exp [fx — fam]2], x— 00. 
Поэтому, умножая обе части равенства (18) Ha x2 и пере- 
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ходя к пределу при х—+ со, получаем 
с, (1) = (1/2md)"e'*4, с, (A) =0. 


Если argh==1, то следует использовать несколько иную 
асимитотическую формулу для Н®; однако окончательный 
вывод остается тем же самым. Мы получаем, таким образом, 
асимптотическое равенство 

ty insta se afin coe need (и) 

НО (ax) == 2х) - зе мвд [1-0 (3), (19) 
где 1ах«ю и ш)>=0. Аналогично доказывается, 
что 


H® (ix) = (1/20) -"зебду, (x) [1 +0 ( =) if (20) 


где 1<Casxtx<oo uw 1) 20. Считая ) положительным и 
полагая х=зесВ, мы получаем после подстановки выраже- 


ний (6) и (14) формулу (2). 


4.5. Точки перехода 


Рассмотрим теперь снова дифференциальное уравнение 


Лиувилля 
Ур РГО y= 0, (1) 


где ) — большое положительное число. Как и в п. 4.1, будем 
считать, что х— вещественное переменное, p(x) — вещест- 
венная функция, имеющая непрерывную вторую производную, 
а г(х) — функция, непрерывная на интервале а=х=6. 
Однако, в отличие от п. 4.1, допустим, что функция р (x) 
может обращаться в нуль в некоторых точках интервала (а, 6). 
Назовем нули функции р(х) точками перехода дифферен- 
циального уравнения (1). Ради определенности предположим, 
что p(x) имеет простой нуль в точке х ==с и не обращается 
в нуль в остальных точках интервала а х=6, причем 
р’ (с) > 0, так что р(х)<0, если axx<c ир(х>0, 
если с«%х=. 

Мы видели выше, что на интервале с-- = х==р, > 0, 
где функция р(х) положительна, решения уравнения (1) 
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асимптотически имеют вид 
с, [р (х)| `"* cos {i \ [p (x)]'2 dx} + 
Не, [p(x)}7"* sin {| К [2 (x)]"? ах} e 2 


Аналогично можно показать, что на интервале а = х=с—в, 
иде функция p(x) отрицательна, решения уравнения (1) асимп- 
тотически имеют вид 


с, [— р (х)| 7" exp {0 { [-— 2) 4х dy 
и as { —1) [-- род] ах] . (3) 


Справедливость этих асимптотических формул связана с тем, 
что функция 2 (2) в формуле 4.1 (4) ограничена, следовательно, 
ни одна из этих асимптотических формул не имеет места 
при х==с. Как показывает равенство (2), справа от точки 
х==с каждое решение уравнения (1) быстро колеблется; 
равенство же (3) показывает, что слева от точки х=с все 
решения уравнения (1) медленно изменяются. Около точки 
X==C происходит переход от одного типа поведения к 
другому. 

Здесь возникают две проблемы. Первая из них состоит 
в «склеивании» асимптотических разложений, то есть в уста- 
новлении связей между постоянными C,, с, и постоянными С,, 
с Которые имеют место, если формулы (2) и (3) являются 
асимптотическими представлениями на различных интервалах 
одного и того же решения уравнения (1); вторая проблема 
состоит в определении асимптотического вида решений урав- 
нения (1) на интервале (с —е, с--®). 

Существуют два остроумных метода решения первой проб- 
лемы. Первый из этих методов был применен Джёффризом 
(1923) и открыт вновь Крамерсом несколькими годами позже. 
Он основан на том, что в окрестности точки х==с функцию 
P(X) можно приближенно заменить линейной функцией 
(х —с)р’ (с), а функцией r(x) можно пренебречь. Получаю- 
щееся при этом дифференциальное уравнение допускает ре- 
шение в бесселевых функциях порядка -- 1/3. Сравнение 
асимптотических формул для бесселевых функций с форму- 
лами (2) и (3) и приводит к искомым формулам, связываю- 


щим коэффициенты с,, с, с коэффициентами с, с.. Второй 
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метод, разработанный Цвааном (1929), заключается в обходе 
точек перехода. Если p(x), г(х) — аналитические функции 
от х, дифференциальное уравнение интегрируется вдоль пути, 
лежащего в комплексной плоскости, который состоит из ве- 
щественных интервалов (a, с—5), (с-|-е, 6) и полукруга 
в комилексной плоскости, соединяющего c—e ис- =. Вдоль 
этого пути функция 0, определенная формулой 4.1 (4), огра- 
ничена, а потому можно применить метод Лиувилля (или 
какой-нибудь вариант этого метода). Этот метод приводит 
к тем же самым формулам склеивания, что и первый. Оба 
метода можно распространить на случай, когда функция 
p(x) имеет нуль произвольного порядка. Они известны как 
WKB- (или иногла WKBS-) методы. Более подробное изло- 
жение этих методов см., например, в книге Морс и Фешбах 
(1960, стр. 90—104). 

Значительно более трудной является вторая проблема, 
а именно, определение асимптотического вида решений урав- 
нения (1) в интервале (с —в, ce). По-видимому, не суще- 
ствует простых элементарных функций, описывающих переход 
от медленных колебаний к быстрым, поэтому приходится 
использовать для получения асимптотических формул неко- 
торые высшие трансцендентные функции. Простейшим диф- 
ференциальным уравнением вида (1), имеющим точку перехода, 
является уравнение 


ЧН Юху ==. (4) 


Решения этого уравнения хорошо известны, и целесообразно 
искать асимитотические формулы для решений уравнения (1), 
выраженные через решения уравнения (4). Преобразование 
Лнувилля 4.1 (2) переводит уравнение 4.1 (1) в дифференци- 
альное уравнение, коэффициенты которого приблизительно 
постоянны. Подобно этому, можно искать преобразование 


E==9(x), п=0у, (5) 


переводящее уравнение (1) в уравнение, приближенно имею- 
щее вид (4). Преобразование (5) переводит уравнение (1) в 


4% 1/9” ф \\ dn рег $ 44, 
le at at [a pea * 


Для того чтобы это дифференциальное уравнение прибли- 


Г 


112 диФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С БОЛЬШИМ ПАРАМЕТРОМ [ГЛ. IV 


женно совпадало с уравнением (4), определим сначала % так, 
чтобы выполнялось равенство 


2 =0. (6) 


Из этого уравнения находим, что Y= 
ляем % из уравнения 


Потом опреде- 


=. (7) 


Если 7) и $ определены указанным способом, то полученное 
выше дифференциальное уравнение принимает вид 


a Ур т, (8) 
где 
19” 3 9" 
p= т — дея. (9) 


Из предположений о функциях р и г, сделанных в начале 
этого пункта, вытекаст, что существует единственная веще- 
ственная функция 9, имеющая непрерывную производную 
третьего порядка и удовлетворяющая уравнению (7). Для 
этой функции точная нижняя грань модуля ее производной 
отлична от нуля, а потому @(§) является ограниченной функ- 
цией. Поэтому можно ожидать, что асимптотическая формула 
для решения уравнения (8) имеет вид 


с, Н, (6) +e, (5), (10) 


где Н, {х) и Н, (х) — линейно независимые решения уравне- 
ния (4). 

Это обобщение метода Лиувилля было виервые применено 
для получения асимптотических формул, выражающих реше- 
ние уравнения (1) на интервале (с—, с-- =); однако ясно, 
что этим методом можно воспользоваться и на всем интер- 
вале (а, 5). При этом вместо трех различных асимптотиче- 
ских формул, имеющих место на интервалах (а, с —®), 
(с—вс--®), (с-[е, 6) соответственно, получается одно 
равномерное асимптотическое представление решений 
уравнения (1) на всем интервале а=х<=6. Этот метод был 
впервые введен Лангером, который разработал его в ряде 
работ, частично указанных в конце этой главы. Среди раз- 
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рабатывавших метод Лангера следует отметить Черри. Cy- 
ществует указатель литературы, касающийся этого метода 


{см. ссылку в конце главы). 
Прежде чем детально описывать этот метод, укажем 
кратко некоторые свойства решений уравнения (4). 


4.6. Функции Эйри 


Дифференциальное уравнение 
dw 
ll eee ee 
ae — 2% 0 (1) 
можно привести к дифференциальному уравнению, которому 
удовлетворяют бесселевы функции порядка 1/3 (Ватсон, 
1949, п. 6.4). Мы будем пользоваться следующими обозначе- 
НИЯМИ;: 
2 asa) : 
c= 2", @ == е?=йз, (2) 
Линейно независимыми решениями уравнения (1) являются 
так называемые Функции Эйри первого и второго вида 
Е 1 ji а = 1 Ид \12 ы 
Ai (2) == 24? [Ty (0) — lus (2)] See (+) Куз <), 
(3) 


2 \12 


Bi) = Pi 9 Ель ©. 
Непосредственное вычисление приводит к равенству 
Ai (2) Bi’ (2) — АР (z) Bi(z)==07 (4) 
Справедливы также формулы 


В 1 ° 
Ai (—2) => z? [Jaa (8) + Jus I, 
m z\t'2 (5) 
Bi(—2)=($)"" Fas) — Jun OL 
При вещественных значениях х имеют место интегральные 
представления 


cos (G+ ) dt, 
(6) 
[exe (—G +t) + sin(F +24) dt. 


114 ДИФФЕРЕНЦИАЛЬНЫЕ УРАВНЕНИЯ С БОЛЬШИМ ПАРАМЕТРОМ [ГЛ. IV 


Эти интегральные представления могут быть преобразованы 
в интегральные представления, которые содержат контурные 
интегралы и справедливы при комплексных значениях х. 
АЁ (2) и Bi(z) являются целыми функциями от 2; при веще- 
ственных значениях 2 они принимают вещественные значения. 
Функция Ai(z) имеет бесконечно много нулей, причем все 
нули этой функции вещественны и отрицательны (Ватсон, 
1949, п. 15.7). Для любого целого m функция w,, (2) = Ai (wz) 
также является решением уравнения (1). Непосредственное 
вычисление показывает, что 


Wry, (2) Wy (2) — W', (2) Wy (2) = 


Яр ("+ +3) | (7) 


Отсюда вытекает, что функции W,, и W,,, линейно незави- 
симы. С другой стороны, любые три w,, из функций связаны 
линейной зависимостью. В частности, 


Wy (2) OW m4, (2) и, , (2) = 0. (8) 
Далее, имеем 
ВЕ (2) = {|9*АЦо*2) — wAi (w2)]- (9) 
Разложения функций Ai(z) и (2) в степенные ряды 
вытекают из равенства (3). В частности, 


3:3 Ai (0) == 3% Bi (0) =—! 


Г (2/3) 
зай ври 1 
pas: =3~-1% = а 
313 Ai’ (0) = 3 ~ 8 ВИ (0) МЕ 

Асимптотическое поведение функции А?(2) было изучено 
ви. 2.6 при —п/3 Cargz<n/3. Полученные результаты 
можно распространить на большую область с помощью пово- 
рота в #-плоскости пути интегрирования в формуле 2.6(3). 
При этом получается 


Аа yn tz e-E[L OK], (11) 


2—0, 9 — a1 <Cargz<n. 


(10) 


Результаты, имеющие место в секторах, содержащих отри- 
цательную полуось, могут быть получены с помощью фор- 
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мулы (8). Они имеют следующий вид: 
1 


Ai (2) == по ее 1 -- 057") -- 1 [1 --067]Ь, 
zo, п/3 < arg z< 51/3; 
Аа тие - 0-9] № [LO 


2— <, — 51/3 Cargz< —n/3. 


(12) 


Из равенства (9) вытекает 
Bi(z) == п-т ле 1-0 (57, (13) 


Zo, —п/3 ап. 


Все асимитотические разложения выполняются равномерно 
по arg 2, если 2 изменяется в замкнутых секторах, целиком 
расположенных внутри указанных выше открытых секторов. 

Из полученных асимптотических формул вытекает, что 


(1 |=“) eS AZ (2), —T<argz<n, (14) 


(l-+z|“)e" Bi(z), —п/3 < ав 2 = п (15) 


являются ограниченными функциями от 2. Ограничены и ве- 
личины, обратные этим функциям, если, разумеется, нули 
функции (14) соответствующим образом исключены; напри- 
мер, это имеет место, если — п-- = ава и— в, => 0; 
функция (15) не имеет нулей. 

Относительно дальнейших сведений о функциях Эйри см. 
книгу Миллера (1946); в ней приведены также таблицы значений 
этих функций. 

Мы закончим этот пункт доказательством следующего 
неравенства. Если wz) является решением уравнения (1), 
‘Удовлетворяющим условиям  (#) =0, w' (1) —1, то при 


р 
x>0, #>0, t=Zx">0, c= 28:0 
имеет место неравенство 


(16) 
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Для доказательства используем теорему сравнения Штурма. 
Функция 
eb-t =" 


20 


9 (х) = 


удовлетворяет дифференциальному уравнению 


“—(х-Нвн)=0 (17) 


и тем же самым начальным условиям в точке в что и функ- 
ция w. Из уравнений (1) и (17) вытекает, что обе функции 
v(x) и W(x) отличны от нуля, если х > 0, #> 0, х=ЕЁ. Рас- 
смотрим функцию f(x) == (х);о (х), х5ЕЁ Очевидно, что 
/(х)—-1, когда x—t, а потому мы положим Л(В=1. 
Функции Vu ® имеют одинаковые знаки при х 522. Поэтому 

а 

ах 


” ” Suw 
(w'v — wo') = w"'v — wo" = oe <0 


и, следовательно, w’U— 0’ является убывающей функцией 
от x. Так как эта функция обращается в нуль при х=, 
то f’ (<) < 0 при x >t uf’ (x) > 0, если Ox <1. Таким 
образом, /” (x) имеет максимум при х = и, значит, Of (x)<1, 
если O<x, Ё< с. Принимая во внимание определение 
функций Л и 9, мы приходим к неравенству (16). 


4.7. Асимптотические решения, справедливые 
в области перехода 


Предположим, что (а, 6) — ограниченный интервал, с — 
внутренняя точка этого интервала, р(х) — вещественная 
функция, заданная на этом интервале, имеющая непрерывную 
производную второго порядка, и такая, что р(х)<0 при 
a<x<c, p(x)>0 при c<x<b, p(c)=0, p'(c) 0. 
Пусть, далее, r(x, 4) — функция, заданная в области а xX <b, 
LES, где $ задается неравенствами || >), 9, Sargh<y,, 
причем функция r(x, \) ограничена в этой области и при 
каждом фиксированном значении A из $ r(x, }) является не- 
прерывной функцией от x на интервале а = х <b. Рассмот- 
рим дифференциальное уравнение 


УНР (x) +r (x, у =0. (1) 
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Обозначим через (x) (единственное) вещественное ре- 
шение уравнения 4.5(7), имеющее непрерывную производную; 
это решение дается формулами 

x 
2 р r = 
Вере (а) == (рб хе, 
с 


(2) 


роду = [[--р 4, exe 


Здесь 9(х)>0, если хе, 9(х)< 0, если x<c H все 
дробные степени имеют положительные значения. 
Функции 
Уи (5%) == [9' (х)]- 2 А -— ет о (x)|, т==0, 1, 


У, (x) = [p" |7 Bi [— "9 (x)] 


zl 


(3) 


удовлетворяют дифференциальному уравнению 
ид [Mp №, 4] у=0, (4) 
где 
” nya 
(Е) 
® $ 
— производная Шварца от $. Это дифференциальное урав- 
нение вытекает из формул 4.6(1) и 4.517). Так как диффе- 
ренциальные уравнения (1) и (4) отличаются сравнительно 
малыми членами, можно рассматривать решения уравнения (4) 
как главные члены формальных решений. 

При сделанных выше предположениях дифференциальное 
‘уравнение (1) обладает решениями, которые в соответ- 
ствующих секторах комплексной \-плоскости асимптоти- 
чески представимы функциями У„ m==—1, 0, 1,2. Дока- 
зательство этого утверждения аналогично доказательствам, 
проведенным в пп. 3.3, 3.4 и 4.3; мы разобьем его на не- 
сколько этапов. 

Пусть У(х) — любое решение уравнения (4), a K(x, t)— 
решение уравнения (4), удовлетворяющее начальным усло- 
вия 


K(t, t)=0, Sa, j=; (5) 
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пусть, далее, Aa<xX,<b. Тогда решение интегрального 
‘уравнения 


хде | Кое AFE DM, y Wat, 6) 
где | 
Е, N= aly их, 0) (7) 


‘удовлетворяет уравнению (1). 

Доказательство. Уравнение (6) является интеграль- 
ным уравиением Вольтерра с непрерывным ядром. Для любого 
фиксированного значения A, принадлежащего области $, суще- 
ствование и единственность решения этого уравнения выте- 
кают из общей теории интегральных уравнений Вольтерра 
(либо могут быть доказаны аналогично тому, как это было 
сделано в п. 4.3). Это решение имеет непрерывную произ- 
водную второго порядка, причем подстановка в уравнение (1) 
показывает, что у удовлетворяет этому уравнению. 

При а=х, t<b для всех \ из 5 выполняется неравен- 
ство 


[Kin NF |= 
Acxp |= ве | — обор" — | ор" | 


Пао we Ol 


» (8) 


где A не зависит om x, t uh. Все дробные степени имеют 
здесь главные значения. 

Доказательство этого неравенства проводится различными 
способами в зависимости от взаимного положения точек x, # 
и в зависимости от arg}. Мы проведем детальное рассужде- 


ние при a<x, #=с, — 1/2 < а 1 == п/2. Так как, в силу 
(3) и 4.6(4), 
У У, — УИ = — 8? (АРВЕ— AiBi') == 18, 
то функция 
K(x, =m * [У, (x) У, (В — У, 0 У, (x)] (9) 


удовлетворяет равенствам (4) и (5). Ho 9<0 и 
[аге ( — №34) | <= т/З, а потому выражения 4.6 (14) и 4.6 (15) 
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ограничены; следовательно, найдется число В, не зависящее 
от x, и», такое, что 


2 
Bexp, — 3 Rel Роз} 


пе в 


[|= 


Аналогичная оценка имеет место и для У,. Пользуясь этими 
оценками и равенством (9), мы убеждаемся в справедливости 
неравенства (8) в рассматриваемом случае. Для того чтобы 
доказать неравенство (8) в остальных случаях, надо выразить 
функции У, и Y,, а следовательно, и функцию А, через два 
линейно независимых решения уравнения (4), одно из которых 
ограничено при k— co (в рассмотренном выше случае таким 
было решение У,) и применить потом оценки, вытекающие 


из 4.6(14), (15). 


Докажем теперь, что существуют решения у„, и =—— 1, 
0, 1, 2, соответствующие решениям Y,, и асимптотически 


представимые ими. При этом оказывается необходимым нала- 
гать в каждом случае ограничения на argh. Эти ограничения 
соответствуют проведению разрезов на \*-плоскости. 

Для Yo и у, мы предположим Ке). 0 и определим эти 
функции как решения интегральных уравнений Вольтерра 


= rir +1 Kes A NPL, Wy, (dt, (10) 


у, (x) = ¥,(x)+- | Ка, NF (t Dy, (Mat. (11) 


Seeks 


Для любого фиксированного значения \ из 5 существование 
решений этих интегральных уравнений вытекает из общей 
теории интегральных уравнений Вольтерра (либо может быть 
доказано методом последовательных приближений). Мы пока- 
жем сейчас, что у, — У» у, —И,, если \— сю, Reh SO. 

Если Ве. =0 и шА=-0 и а=х=ь или если } >=0 
и ах = с, то У, =5=0, У, 20. В этом случае можно поло- 
жить у, = 7,2, у, ==7,2, и получить интегральное уравнение 
для функций 2, и 2,. Так как при сделанных предположениях 
функции 4.6(14) и 4.6(15) ограничены и точные нижние грани 
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их модулей отличны от нуля, то 


LAG) ety] LIT" + Lg tl" 
Yee [Е a | ее 


X exp ( — 3-е {[ Wp (H]" — [25969]. 


Аналогичная оценка имеет место для | Y, (1) /У, (х) |. Комбинируя 
эти оценки с неравенством (8), мы получаем неравенства, 
которым удовлетворяют ядра интегральных уравнений для 
2, и 2,. Они имеют вид 


С 


У, ,. |- ee: . 

ak. ВР, | «тит (12) 
axt<x<bu шт) 520 или аз хжси > 0. 

У, (t) ‘ С 

Pak OF | < ртути! (13) 


ax<x<b, |t—c|<|x—c| и ш\520 
или ах си \ > 0. 


Для каждого фиксированного ) интегральное уравнение для 
2„ m==0,2, имеет ограниченное решение. Пусть Z,,(d) 
является максимумом для |2„(х)| при a<x<b. Из интег- 
ральных уравнений и неравенств (12), (13) имеем 


b 
о a at 
Z, 0) <1+4+22Z, wl ae 


Для достаточно больших значений |)| выполняется неравен- 
ство |Z,,())|<<2, а из соотношений (10), (11), (12), (13) 
вытекает, что если Ве). =0, а=х= В, и при вещественных 
значениях h имеем @<x<c, то 


Im (X) = Yip (x) [1- 007] (14) 


равномерно по х, когда h—+ со, m==0,2. 

В случае, когда \ принадлежит положительной полуоси, 
функции У, и У, имеют нули при x >> ¢ и соотношение (14) 
теряет силу вблизи этих нулей. Применимый в этом случае 
результат может быть получен из соотношений (10) и (11) 
< помощью равенства (9) и оценок для У, (А) и Y,(¢). Из 
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равенства (9) получаем 
y= У, (ху А (У, (ВР, Dy, да — 


—mh~ "HY, (x) [YF Ny, Oat 


a 


Разобьем эти интегралы следующим образом: 


ee) 


и применим к первому интегралу оценку (14). Мы получим 
Yo (=, (x) [1 + OO" +- Y, (x) 97") (15) 


равномерно по x, когда h—+ oo, cM X <b>. 
Точно Tak же доказывается, что 


У: ==, (x) [1 + OG") + У, ©) O07) (16) 


равномерно по x, когда h—>+ oo, cx x <b, \>0. Равенства 
(14) — (16) полностью описывают асимптотическое поведение 
решений у, и у, 

Совершенно аналогично исследуются решения ys. Мы 
определяем у, иу_, с помощью интегральных уравнений 


ь 
In )=Vq (x) —) K(x, DF у, (04, m=1, —1 (17) 


и исследуем решения этих уравнений точно так же, как 
исследовали решения уравнений (10), (11); при этом в случае 
решений у, предполагаем 1-0, а в случае решения у_, 
предполагаем Пи). = 0. Функции У, и Y_, имеют нули, если 
\ — чисто мнимое число и х< с. В этом случае асимптоти- 
ческие формулы надо несколько видоизменить. Окончатель- 
ные результаты имеют следующий вид: 


у. =И, (x) [1-007] (18) 


равномерно по xX, когда h—->o, axx<b, ImiSO, 
Ве} 520 либо ex x<b, —A>0. 


AM=YC)M+OAM+Y.@)OR) (19) 
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равномерно по х, когда \ — oo, а=х=с, —Йй>0. 
y,W=Y_,@%)[1+007) (20) 


равномерно по х, когда \ —+ со, а x <b, Im) <0, ReA~0 
либо с=х=, th>0. 


У (=, ПО] У, (OG) (21) 


равномерно по x, когда ee axx<c, h>O. 

Равенства (14)— (16) + г (18) — (21) и доказывают спра- 
ведливость результата, сформулированного в начале этого 
пункта. С помощью более детального исследования инте- 
гральных уравнений доказывается, что ул — Ум. 

Если функцию r(x, 2) можно разложить по степеням \-', 
то такое разложение возможно и для формальных решений 
уравнения (1). Приближения, полученные в этом пункте, 
можно рассматривать как главные члены формальных реше- 
ний. Как это было и выше, существуют два типа формаль- 
ных решений. Первый из них соответствует типу 4.23) и 
имеет вид 

n 


¥ (x) Ха," + (x) DB (HRM, (22) 


n 


где У— решение уравнения (4), a a, (x) и B, (x) — функции 
от x, не зависящие от hk. Рекуррентные дифференциальные 
уравнения для этих функций получаются путем подстановки 
разложения (22) в уравнение (1) и сравнения коэффициентов 
при одинаковых степенях ).. Это приближение было применено 
Лангером (1949). 

Второй тип формальных решений соответствует решениям 
4.29). Эти решения имеют вид (3), за исключением того, 
что вместо функции ф(х) мы имеем функцию 9 (x, 4), которая 
зависит от A и обладает формальным разложением 

[2] 


Хы" (23) 


n=0 


В этом разложении ©, (x) — дифференцируемое вещественное 
решение уравнения 4.5(7), а функции %,, $„, ... удовлетво- 
ряют рекуррентным дифференциальным уравнениям, получае- 


4.8] ПРЕДСТАВЛЕНИЯ БЕССЕЛЕВЫХ ФУНКЦИЙ 123 


мым путем подстановки функции 
gy thAi(— Wg) 


в уравнение (1) и сравнения коэффициентов при одинаковых 
степенях hk. Такие решения были использованы Черри (1950). 
Дифференциальное уравнение 


Уч DY=0, Ge N= Day ye 


является более общим, чем уравнение (1), поскольку оно со- 
держит член 4, (х)}. Из-за этого члена возникают дополни- 
тельные осложнения, которые были исследованы Лангером 
(1949). 

Обобщение этих результатов на случай комплексного пере- 
менного было предпринято Лангером (1932) и Черри (1950); 
случай, когда переменное x изменяется в бесконечном интер- 
Bane, был исследован Черри (1950). 


4.8. Равномерные асимптотические представления 
бесселевых функций 


В заключение этой главы применим результаты предыду- 
wero пункта к дифференциальному уравнению 4.4{3). С ux 
помощью мы получим асимптотическое представление для 
Л (ix), справедливое при \—+ сю, Reh sO равномерно для 
всех положительных х. Эти результаты содержат результаты 
п. 4.4, относящиеся к функции Л (kX) и, кроме того, запол- 
няют пробел, оставшийся в п. 4.4, а именно случай р ха 
(< 1<а. 

Итак, применим методы предыдущего пункта к уравнению 
4.4(3), для этого уравнения 


1 
рот, r(x, N= a) 


и потому точкой перехода является х==1. Функция о опре- 
деляется дифференциальным уравнением 


(2) 
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и из равенства 4.7(2) вытекает, что 


О fet 1y'sat= 


=, Oe, (8) 


Зо: =[u- 17") 341 = 
== (х° — 1)'2 —arecosx7'= (хх), 1<x<oo (4) 


(см. также 4.4(4) и (14)). 
Так как G(x) является аналитической функцией и 
@' (x) 5 0, = очевидно, функция 


Lol” 3 fo" \? 1 
Роду ors 1 (1) 4 


непрерывна на интервале 0 х< с. Для того чтобы изу- 
чить поведение функции F(x) при х—+0 и х—+ oo исполь- 
зуем формулу для производной Шварца от сложной функции. 
Она имеет вид 


te, ха и (Z) tu >}. 


Положим в этой формуле и=3 при х<1 и u==f при 
x > 1. Простое вычисление показывает, что 
53” [8 аа 
P(x) = Tor +5 ae = Ter Tia 
0 ‘ ХВ, 


5 4-х 
РЕВ 1<я<. 


Из (3) и (4) вытекает 
B==O(Inx), ==0(х7'), x—0, 
1=0(х), f= O(1), x0, 
а потому 
F (x)= O[(xInx)7*],. x0, 


F(x) =O(x7*), х— о. (5) 
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Кроме того, как и в п. 4.7, получаем неравенство 


У, (0 С 
Dlg Olewrs =a la (6) 


теряющее силу вблизи нулей функции V(x). 

Мы видим, таким образом, что точку х==0 можно при- 
нять в качестве постоянного предела в нитегральном урав- 
нении для 2, (x). Из соотношений (5) и (6) вытекают оценки 
для ядра этого уравнения, которые показывают, что это ядро 
ичтегрируемо на интервале (0, со). Несмотря на то, что 
x==0 является особой точкой дифференциального уравнения, 
можно применить метод п. 4.7. При т ==0 справедливо урав- 
нение 4.7(14); при этом, так как из соотношений (5) и (6) 
i in. зто 


x 
IY, (2) ЧЕ = 1 
(xe oe оо ( nap) = (mz): а 
© 


то при малых значениях x член О(\`') 4.7(14) можно заме- 
нить более точным выражением 


O( rm): (7) 


Мы получаем, таким образом, что у дифференциального урав- 
нения 4.4(3) существует решение у, (x), для которого равенство 


У, (x) =’ 7 2 АЦ — $) fl + 00-1) (8) 


выполняется при 0% хх, А+, Ве’ 0, за исключе- 
нием случая, когда xX > Ги } — положительное вещественное 
число; в этом случае надо несколько изменить вид остаточ- 
ного члена. Пр OMe <b< 1 остаточный член можно за- 
менить более точным выражением (7). 

Для того чтобы выразить у, через бесселевы функции, 
положим в равенстве (8) 0 < x <b < 1; тогда $ (х) = 9 (5) < 0, 
функцию Эйри можно заменить асимптотическим представле- 
нием 4.6(11), О-член можно заменить более точным выраже- 
нием (7), и мы получим 


о) = TT 8 (— уе [1-50 (==) | = 
п- 1) - “eae? +0(т=)|, 0<x<b<l, 


K(x, Е |= 


ы|- ы- 
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где a — функция, определенная формулой 4.4(5). Сравнивая 
это выражение с 4.4(6) и 4.4(13), получаем 
Qe ки а 
ле Ra у, 0х. 
Следовательно, 


t 


2h 


J. Qe) = EHO PT" eg!) "2 AL(— Hg) [1+ OO]. 9) 


Для этого равенства справедливы те же замечания, касаю- 
щиеся остаточного члена, что и в равенстве (8). Применяя 
к TQ-+1) формулу Стирлинга, можно представить этот 
результат в упрощенной, но ослабленной форме 


Ji, (1x) = (ях) АНИ -О(-*)]. (10) 


Это равенство имеет место равномерно по x, O<x< oo, 
когда \-—* сою, Reh 0, причем остаточный член надо не- 
сколько изменить вблизи нулей функции Al ( — 2s). Заметим, 
что в формуле (10) остаточный член содержит остаточный 
член формулы Стирлинга и не может быть уточнен при ма- 
лых значениях х. 

Полученный результат содержит в качестве частного слу- 
чая формулу 4.4(Т). Покажем, что равенство (19) содержит 
в качестве частного случая и формулу, получаемую сумми- 
рованием обеих формул 4.4(2). Предположим, что x >= a>), 
g(x) >= 9 (а) > 0 и применим 4.6(12). Мы получаем 


J, 0.x) ~ (+ nix) Fy! = tae ~ "he cos [Af (x) — 0/4], 
1«а<х«о, 
или 
г . \7 ‚ 
Jy (Ax) ~ (zm) a(x) cos [if (x) — 1/4], 
м / 
1«а=х«о, 


что согласуется с формулами 4.4(2). 

Основной результат этого пункта (9) был распространен 
на комплексные значения х Черри (1948), который получил 
также приближения высшего порядка, 
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